AD-A072  094  GEORGETOWN  UNIV  WASHINGTON  0 C DEPT  OP  PHYSICS  F/6  20/1 

INFLUENCE  OF  ABSORPTION  ON  NON- SPECULAR  ULTRASONIC  REFLECTIVITY— ETC (U) 
JUL  79  TO  N60C*  W G WATER  N00014-78-C-0584 

UNCLASSIFIED  GUUS-07791  NL 


DDC  FILE  COPY  ABA072094 


SECURITY  CLASSIFICATION  QF  This  RAGE  rMTiwi  Data  Knl.t.rf) 


READ  INSTRUCTIONS 


REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 


S CATALOO  NUMBER 


[I.  OOVT  ACCESSION  NO 


N0t?Q14-78  -C-0584]f 


Physics  Department ' t 
Georgetown  University 
Washington  DC  20057 


121108 


CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research, 
Arlington,  Virginia 


IS.  SECURITY  CLASS,  ft  tftU  report; 


MOWlTORl NO  AGENCY  NAMt  i AQPRlftVJI  different  free i Cewtrelllnj  Off! cej 


Unclassified 


SCHEDULE 


Ift.  blSTRltUtlON  STAtCMCNT  (ml  (fill  ftepert) 


Approved  for  public  release;  distribution  unlimited 


Approved  for  public  release;  distribution  unlimited 


It.  REV  WOROS  ICanllnua  an  rararaa  ilFa  It  naaa.aary  ana  INanllfF  *F  Man*  numSa*) 


ultrasonic  reflectivity,  bounded  beam,  numerical  integration 
method,  longitudinal  critical  angle,  Rayleigh  critical  angle, 
ultrasonic  absorption 


A numerical  integration  method  is  developed  for  calculating  the 
changes  in  beam  profile  and  amplitude  upon  reflection  of  a boun- 
ded ultrasonic  beam  Incident  on  a liquid-solid  flat  interface  at 
the  longitudinal  wave  critical  angle  and  the  Rayleigh  angle.  No 
approximations,  valid  only  at  the  Rayleigh  angle,  are  used,  and 
the  method  is  applicable  to  all  angles  of  incidence. 


Unc lassified 


This  Technical  Report  describes  work  performed  by 
T.  D.  K.  Ngoc  for  his  doctorate  in  Physics  at  George- 
town University.  The  study  contained  in  the  Report 
is  part  of  the  Project  supported  by  the  Office  of 
Naval  Research,  Physics  Program  (Code  421),  under 
Contract  N00014-78-C-0584. 

The  study  is  concerned  with  the  influence  of  absorp- 
tion on  the  numerical  values  of  the  reflectivity  of 
a bounded  ultrasonic  beam  reflected  from  a liquid- 
solid  interface.  Particular  attention  is  given  to 
changes  in  the  beam  profile  upon  reflection  at  the 
critical  angles.  Only  approximate  values  of  the 
beam  parameters  had  been  available  until  now;  the 
numerical  integration  method  described  in  this  Re- 
port yields  exact  values  for  all  angles  of  incidence. 
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ABSTRACT 

A theoretical  and  experimental  investigation  of  non-specular 
ultrasonic  reflectivity  at  a liquid-solid  interface  is  described. 

The  role  of  attenuation  in  the  formulation  of  the  plane  wave  reflection 
coefficient  is  analyzed.  This  attenuation-dependent  reflection  co- 
efficient is  used  in  a numerical  integration  method  which  is  developed 
for  the  calculation  of  the  reflected  beam  profile.  The  results  are 
compared  with  other  existing  formulations.  The  reflected  profile  Is 
calculated  by  this  numerical  approach  for  two  specific  cases:  inci- 
dence near  the  Rayleigh  angle  for  a water-stainless  steel  interface 
and  incidence  near  the  longitudinal  critical  angle  for  a water-Plexl- 
glas  interface.  In  each  case,  the  non-specular  reflection  effects 
are  investigated  as  one  parameter  is  varied  and  the  others  are  kept 
unchanged.  Comparisons  are  made  with  results  obtained  by  the  Bertonl 
and  Tamlr  formulation  for  incidence  at  the  Rayleigh  angle.  The  theo- 
retical predictions  are  then  compared  to  experimental  results  obtained 
by  the  Schlieren  visualization  technique. 
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CHAPTER  I 


INTRODUCTION 

The  reflection  coefficient  derived  from  the  plane  wave  theory  is 
inadequate  in  describing  the  reflection  of  ultrasonic  waves  Incident 
from  a liquid  medium  onto  a flat  liquid-solid  interface.  In  practice, 
the  sound  waves  are  normally  generated  by  a source  of  finite  dimensions 
and  thus  cannot  be  considered  as  plane  waves.  From  the  plane  wave 
theory,  it  is  expected  that  the  reflected  sound  beam  has  the  same  in- 
tensity distribution  as  that  of  the  incident  beam  at  the  reflecting 
Interface;  or,  in  other  words,  it  is  specularly  reflected.  However, 
several  interesting  phenomena,  known  as  non-specular  reflection  effects, 
have  been  observed  in  experiments  as  a result  of  the  boundedness  of  the 
incident  beam,  especially  for  incidence  at  the  Rayleigh  critical  angle. 

The  Raleigh  critical  angle  is  a particular  angle  of  incidence  at 
which  the  incident  longitudinal  wave  is  mode-converted  into  a leaky 
Rayleigh  surface  wave.  The  characteristics  of  thts  surface  wave  are 
that  its  magnitude  decreases  rapidly  with  distance  from  the  interface. 
Its  propagation  velocity  is  less  than  the  shear  wave  velocity  in  the 
solid,  and  it  loses  energy  through  reradiation  into  the  liquid.  For  a 
liquid-solid  structure,  the  plane  wave  theory  also  predicts  the  exis- 
tence of  the  longitudinal  and  shear  critical  angles.  The  longitudinal 
critical  angle  is  the  angle  of  incidence  beyond  which  only  refraction 
of  shear  waves  is  possible;  and  the  shear  critical  angle  is  the  lower 
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limit  of  th«  incident  angular  range  that  results  In  total  reflection. 

The  first  analytical  formulation  developed  to  treat  the  reflec- 
tion of  a bounded  ultrasonic  beam  for  a liquid-solid  interface  was  es- 
tablished by  Schoch  [l] . Schoch's  model  predicts  a lateral  displace- 
ment of  the  reflected  beam  for  all  Incident  angles,  which,  including 
the  Rayleigh  angle,  result  In  total  reflection.  Nevertheless,  his  model 
fails  to  account  for  other  non-specular  reflection  features,  such  as 
the  possible  intensity  null  and  the  widening  of  the  reflected  beam,  ex- 
perimentally observed  for  Incidence  at  the  Rayleigh  angle.  Another 
formulation,  derived  by  Bertonl  and  Tamir  [2]  to  treat  non-specular  re- 
flectivity specifically  for  Rayleigh  angle  Incidence,  eliminates  these 
discrepancies  and  explains  the  basic  non-specular  reflection  phenomena 
at  that  critical  angle. 

In  essence,  these  two  formulations  are  based  upon  the  principle 
of  spectral  representation  of  a bounded  beam.  Their  main  difference 
stems  from  their  different  methods  of  approximating  the  reflection  co- 
efficient. Schoch  approximates  the  reflection  coefficient  as  a product 
of  its  magnitude  and  phase  with  the  former  being  assumed  constant  and 
the  latter  slowly  varying.  The  expression  of  the  reflection  coefficient 
is  a ratio  whose  numerator  and  denominator  are  complex  functions  of  the 
Incident  angle  and  both  become  aero  for  Incidence  at  the  Rayleigh  crit- 
ical angle.  Taking  advantage  of  the  existence  of  a pole-zero  pair  at 
the  Rayleigh  angle,  Bertoni  and  Tamir  used  a complex  Laurent  series  ex- 
pansion to  simplify  the  mathematically  complicated  expression  of  the 
reflection  coefficient  In  terms  of  its  pole  and  zero.  It  Is  then 
evident  that  the  Bertoni  and  Tamir  formulation  describes  non-specular 
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reflection  phenomena  which  are  caused  by  the  interference  of  the  specu- 
larly reflected  beam  and  a trailing  field,  that  is,  a widening  of  the 
reflected  beam.  Schoch,  on  the  other  hand,  explains  only  a beam  dis- 
placement resulting  from  the  phase  shift  generated  by  the  reflection  of 
sound  waves  at  the  interface.  Because  of  Its  incomplete  assumptions 
concerning  the  reflection  coefficient,  Schoch* s approach  has  not  pro- 
vided correct  results,  especially  at  the  critical  angles. 

The  Bertoni  and  Tamir  formulation,  although  providing  predictions 
qualitatively  consistent  with  experimental  observations,  has  some  limi- 
tations. It  Is  applicable  only  at  the  Rayleigh  critical  angle  and 
leads  to  analytical  difficulties  when  extended  to  non-Gauss ian  profiles 
for  the  incident  beam.  Nonetheless,  its  success  is  one  of  the  reasons 
why  non-specular  reflectivity  at  the  Rayleigh  angle  has  been  studied 
extensively.  Since  the  Bertoni  and  Tamir  theory  applies  only  to  Ray-  - 
lelgh  angle  incidence,  little  is  known  about  the  possible  non-specular 
reflection  phenomena  for  bounded  beams  Incident  at  or  near  the  longi- 
tudinal and  shear  critical  angles.  Preliminary  experimental  evidence 
Indicating  the  existence  of  non-specular  reflectivity  similar  to  that 
at  the  Rayleigh  angle  was  reported  by  Neubauer  [3]  for  Mallory  1000  at 
the  shear  critical  angle  and  for  Pyrex  glass  at  the  longitudinal  crit- 
ical angle.  However,  there  exists  no  published  theoretical  analysis  of 
these  results. 

The  difficulties  in  a theoretical  analysis  of  the  possible  non- 


specular  reflectivity  near  the  longitudinal  and  shear  critical  angles 
are  mainly  due  to  the  singular  behavior  of  the  reflection  coefficient 
at  these  angles.  The  analytical  properties  of  the  reflection 
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coefficient,  concerning  Its  complex  singularities,  together  with  the 
contour  Integration  process,  involved  In  the  spectral  representation  of 
a bounded  beam,  have  been  studied  extensively  by  several  authors  [2,4,5] 
for  the  lossless  situation.  Inclusion  of  attenuation  in  the  media  leads 
to  a more  realistic  description  of  the  non-specular  reflection  effects 
and,  at  the  same  time,  removes  these  analytical  difficulties  as  will  be 
shown  later. 

Attenuation  in  the  media  further  complicates  the  physical  problem 
under  consideration.  Becker  and  Richardson  [4]  investigated  the  in- 
fluence of  absorption  on  reflectivity  to  find  that  the  reflection  co- 
efficient behaves  differently  when  absorption  is  taken  into  account. 
Total  reflection  is  then  not  expected  at  any  angle  of  incidence  and  the 
reflection  amplitude  and  phase  are  shown  to  significantly  deviate  from 
the  lossless  case.  Bertoni  and  Hou  [6]  have  recently  incorporated  at-, 
tanuatlon  into  a theoretical  investigation  of  non-specular  reflectivity 
at  the  Rayleigh  angle.  Pitts,  Fiona,  and  Mayer  [5]  have  extended  Ber- 
toni and  Tamlr's  model  to  a more  complicated  case  of  a solid  plate 
Immersed  In  a liquid.  In  the  same  analytical  framework,  Breazeale  et  al 
[7]  examined  the  reflected  beam  profiles  observed  at  large  distances 
from  the  liquid-solid  interface. 

This  thesis  will  deal  with  the  problem  of  non-specular  reflectivi- 
ty In  the  presence  of  attenuation  covering  the  complete  range  of  inci- 
dent angles.  The  reflected  beam  profile  will  be  calculated  via  a numer- 
ical integration  method  for  a Gaussian  incident  beam  to  facilitate  com- 
parison with  earlier  calculated  results.  Finally,  the  theoretically 
obtained  profile  will  be  experimentally  verified.  A Schlteren 
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visualization  technique  is  used  for  the  experimental  verification. 

Chapter  II  is  a presentation  of  a detailed  investigation  of  the 
plane  wave  reflection  coefficient  for  the  liquid-solid  case.  The  re- 
flection coefficient  plays  an  important  role  In  determining  the  re- 
flected beam  profile  since  the  bounded  incident  beam  is  to  be  considered 
as  a superposition  of  an  infinite  number  of  plane  waves  each  being  re- 
flected In  accordance  with  the  plane  wave  reflection  coefficient.  The 
influence  of  absorption  in  the  non-specular  reflectivity  problem  is 
considered  by  explicitly  Introducing  the  attenuation  parameters  into 
the  reflection  coefficient.  Special  emphasis  is  given  to  the  behavior 
of  the  reflection  amplitude  and  phase  near  the  critical  angles. 

in  Chapter  III,  the  existing  theoretical  methods  of  treating  non- 
specular  reflectivity  are  discussed  in  detail.  A condensed  review  of 
the  Schoch  and  the  Bertoni  and  Tamir  formulations  is  presented  in  prep- 
aration for  the  introduction  of  the  numerical  integration  method.  The 
chapter  is  concluded  with  a comparative  analysis  of  these  formulations. 

Chapter  IV  ts  devoted  to  a presentation  and  discussion  of  the  re- 
sults of  the  numerical  integration  method,  it  was  observed  [2]  that  the 
reflected  beam  profile  is  a strong  function  of  the  angle  of  incidence, 
the  sound  frequency,  and  the  beam  width.  The  most  significant  contri- 
bution of  the  numerical  integration  method  is  its  ability  to  describe 
the  reflected  profile  for  all  angles  of  incidence.  In  addition,  when 
attenuation  is  included,  the  frequency  dependence  of  the  reflected  pro- 
file is  directly  established  in  this  formulation.  The  reflected  beam 
profile  is  then  calculated  for  different  values  of  one  parameter  while 
keeping  the  others  constant.  Although  the  calculation  can  be  performed 
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for  any  angle  of  Incidence,  the  results  given  here  are  limited  to  angular 
ranges  near  the  critical  angles  where  non-specular  reflectivity  Is  prom- 
inent. These  results  are  compared  to  those  obtained  by  the  Bertonl  and 
Tamlr  method  wherever  that  method  Is  applicable. 

Experimental  results  are  presented  and  analyzed  in  Chapter  V. 

Due  to  the  broad  nature  of  the  non-specular  reflectivity  problem,  the 
experimental  effort  Is  focused  only  on  Its  major  aspects.  Variation  of 
the  reflected  beam  profile  as  a function  of  the  angle  of  incidence  near 
the  Rayleigh  and  longitudinal  critical  angles  is  investigated  for  water- 
Plexiglas  and  water-stainless  steel  interfaces. 

The  final  chapter  contains  a summary  of  the  main  results,  possible 
extensions  of  the  numerical  integration  approach,  and  questions  remain- 
ing to  be  answered. 

During  the  course  of  the  theoretical  investigation,  four  computer 
programs  written  In  the  FORTRAN  IV  language  were  developed  to 

(1)  Calculate  the  magnitude  and  the  phase  of  the  plane  wave  re- 
flection coefficient, 

(II)  Search  for  the  locations  of  the  Rayleigh  pole-zero  pair, 

(.tit)  Determine  the  reflected  profile  using  Bertoni  and  Tamlr’s 
results,  and 

(Iv)  Compute  the  reflected  profile  through  the  numerical  inte- 
gration process. 

Representative  versions  of  these  programs,  which  are  at  times 
slightly  modified  to  suit  certain  calculation  purposes,  are  listed  for 
reference  In  the  appendices. 
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PLANE  WAVE  REFLECTION  COEFFICIENT 

i 

FOR  AN  ABSORPTIVE  LIQUID-SOLID  INTERFACE 

* 

The  plane  wave  reflection  coefficient  plays  a central  part  In  de- 
termining the  reflected  profile  for  a bounded  Incident  beam.  Under- 
standing its  behavior  with  respect  to  different  angles  of  incidence 
provides  a useful  perception  of  the  non-specular  reflectivity  problem. 

This  chapter  examines  the  singularities  of  the  complex  plane  wave  re- 
flection coefficient  and  the  variation  of  its  amplitude  and  phase  as 
functions  of  the  angle  of  incidence.  The  first  section  describes  in- 
vestigations of  the  analytical  properties  of  the  lossless  reflection 
coefficient  for  various  typical  liquid-solid  combinations.  In  the  se- 
cond section,  various  published  treatments  of  attenuation  are  reviewed 
and  a specific  manner  of  incorporating  absorption  into  the  reflection 
coefficient  is  presented.  With  attenuation  included,  the  final  section 
is  devoted  to  a detailed  investigation  of  these  aspects  near  the  three 
critical  angles. 

A.  Lossless  Reflection  Coefficient  for  a Liquid-Solid  Interface 

Consider  an  infinite  longitudinal  plane  wave  with  an  angular  fre- 
quency (i),  incident  from  a liquid  medium  at  an  angle  8.,  and  onto  a li- 
quid-solid interface  of  two  isotropic,  perfectly  elastic  infinite  half- 
spaces. The  incident  wave  travels  with  a velocity  v and  part  of  its 
energy  Is  reflected  back  into  the  liquid.  The  other  part  is  refracted 
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as  two  separate  waves  traveling  In  the  solid  with  a longitudinal  wave 
velocity  v^  and  a shear  wave  velocity  vs>  respectively.  The  coordinate 
system  is  chosen  such  that  the  x-axIs  Is  along  the  interface  and  the 
z-axis  normal  to. It.  There  Is  no  y variation  since  the  wave  is  assumed 
to  be  Infinitely  extended  in  the  y direction.  The  above  description  is 
illustrated  in  Fig.  1.  The  derivation  of  the  amplitude  reflection  coef- 
ficient for  such  ;>  liquid-solid  Interface  has  been  carried  out  in  detail 
by  several  authors  [1,8,9].  A convenient  form  of  the  plane  wave  reflec- 
tion coefficient  for  the  lossless  case  is  given  [9]  by 


R(k  ) 
x 


(ks  ' 2kx)2  + kkl*s*d  " pksKd/K 
y y 1 o L * 
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p ■ liquid  density/solid  density, 
k « w/v, 
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The  function  R(k  ) is  complex  due  to  the  existence  of  several 


square-root  terms  In  (1).  In  general,  the  analytical  behavior  of  R(k^) 
is  functionally  determined  by  its  singularities  in  the  complex  k^-plane 
and  hence  should  be  studied  first  of  all  by  considering  these  singular- 
ities. For  the  lossless  case,  It  Is  well  known  that  there  often  exist 


article  displacement  with  single  arrows. 


three  pairs  of  branch  points  on  the  real  nxls  located  at  +kd,+k^,  end 
+k  and  a pole-zero  pair  kp-kQ  positioned  slightly  off  the  real  axis  as 
shown  for  a typical  liquid-solid  interface  In  Fig.  2.  Under  the  boun- 
dary conditions  associated  with  the  lossless  case,  it  can  be  shown  that 

* 

kQ  - kp.  Physically,  these  singular  points  are  related  to  the  critical 
angles  via  the  following  relationships  (Snell's  law) 

smed  - ^ i Sine5  - . smer  - — , (2) 

where  9^,  9#,  and  9^  are  the  longitudinal,  shear,  and  Rayleigh  critical 
angles,  respectively. 

Another  way  to  gain  a better  understanding  of  the  plane  wave  re- 
flection coefficient  Is  to  Investigate  the  amplitude  and  phase  of 
R(kx).  The  behavior  of  the  amplitude  and  phase  strongly  influences  the 
reflection  character Istics  of  the  non-specular  reflectivity  problem  to 
be  considered  later.  A simple  computer  program  (Appendix  A)  was  written 
to  numerically  calculate  the  magnitude  and  phase.  A particular  pro- 
gramming problem  appears  in  dealtng  with  the  square-root  terms  in  the 
expression  for  R ( k^ ) . In  order  to  obtain  a consistent  phase  factor  for 
R(k^),  the  magnitude  and  phase  of  these  terms  are  determined  Individual- 
ly before  they  m«  Insetted  Into  the  calculation  of  the  expression  (I). 

The  calculated  results  are  presented  in  graphical  form  In  Figs. 

3 and  it  for  two  Interfaces:  water-sto Inless  steel  and  water-PlexIglas. 


The  values  of  sound  velocities  and  density  ratios  used  for  these  mater- 
ials era  given  In  Table  I.  The  curves  for  the  water-stainless  steal 
case,  In  general,  represent  a typical  water-metal  combination.  Total 
reflection  at  the  longitudinal  critical  angle  and  beyond  the  shear 
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TABLE  1 


SOUND  VELOCITIES  AND  OENSITY  RATIOS 
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critical  angle  occurs  as  is  usually  predicted  for  lossless  media.  The 
phase  curve  exhibits  a small  dip  at  the  longitudinal  critical  angle  and 
an  almost  360°  phase  shift  at  the  Rayleigh  angle.  Special  variation  one 
might  expect  at  the  shear  critical  angle  is  not  apparent  due  to  its  close 
proximity  with  respect  to  the  Rayleigh  critical  angle. 

Because  the  shear  wave  velocity  of  Plexiglas  is  less  than  the 
sound  velocity  in  water,  the  magnitude  and  phase  of  the  reflection  co- 
efficient  for  the  water-Plexiglas  interface  behave  in  a different  man- 
ner compared  with  most  other  liquid-solid  combinations.  There  exist 
only  two  pairs  of  branch  points  corresponding  to  the  longitudinal  crit- 
ical angle  and  the  angle  of  glancing  incidence,  0j  ■ 90°.  The  Rayleigh 
pole-zero  pair  is  suppressed  and  no  Rayleigh  surface  wave  can  be  gener- 
ated in  this  case.  The  magnitude  of  the  reflection  coefficient  has  a 
very  sharp  peak  at  the  longitudinal  critical  angle  and  decreases  to  al- 
most zero  before  gradually  returning  to. total  reflection  at  glancing 
incidence.  Starting  from  the  longitudinal  critical  angle,  the  phase 
drops  slowly  from  360°  to  180°  at  0j  - 90°.  The  manner  in  which  the 
phase  decreases  in  this  range  is  found  to  depend  mainly  on  the  differ- 
ent values  of  the  shear  wave  velocity  used  for  different  Plexiglas 
samples.  It  should  be  noted  that  the  slopes  of  the  magnitude  and  phase 
curves  exhibit  a discontinuity  whenever  there  exist  a pair  of  branch 
points. 

B.  Incorporation  of  Absorption  into  Reflection  Coefficient 

It  has  been  realized  [2,4,6]  that  attenuation  in  materials  would 
affect  the  non-specuiar  reflection  phenomena  at  the  Rayleigh  critical 
angle.  The  current  procedure  to  account  for  attenuation,  which  has 
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proved  to  be  satisfactory,  Is  to  define  a loss  parameter  from  the  tra- 
ditional attenuation  constant  and  to  Incorporate  It  In  the  wave  number, 

which  then  becomes  complex.  Denoting  these  loss  parameters  as  a,  a. 

o 

•nd  a#  for  the  sound  waves  In  the  liquid,  the  longitudinal  and  shear 
waves  In  the  solid,  respectively,  one  can  write  the  complex  wave  num- 
bers as 

k • (w/v)  (I  + ia/2n)  , 

kd  " (w/vd>  + >«d/2ir)  * (3) 

k#  - (u/vs)  (I  •*>  las/2ir). 

The  loss  parameters  can  then  be  identified  as  attenuation  per 
wavelength  In  nepers  and  are  related,  for  example,  to  the  attenuation 
constant  o of  the  liquid  medium  by 


a 


a\  , 


(*) 


where  f and  X are  the  sound  frequency  and  wavelength. 

In  investigating  the  influence  of  attenuation  on  the  reflection 
coefficient,  earlier  authors  [2,1*]  observed  that  Mk^)  depends  on  the 
loss  parameters,  or  the  product  aX,  rather  than  on  attenuation  or  fre- 
quency alone.  Accordingly,  they  made  use  of  different  values  of  the 
loss  parameters  without  considering  separately  either  frequency  or  at- 
tenuation. However,  some  results  obtained  in  such  calculations  are  of 
little  practical  use  since  attenuation  in  most  materials  is  a non- 
linear function  of  frequency.  Thus  one  does  not  have  the  freedom  to 
separately  select  the  values  of  frequency  and  attenuation  to  arrive  at 


a given  set  of  loss  parameters.  In  practice.  It  is  more  desirable  to 
be  able  to  treat  frequency  as  an  Independent  quantity.  This  can  be 
achieved  provided  the  functional  relationship  between  attenuation  and 
frequency  can  be  determined. 

Another  aspect  concerning  the  method  of  introducing  attenuation 
into  the  reflection  coefficient  by  the  same  authors  [2,4]  should  also  be 
mentioned  here.  In  dealing  with  reflection  phenomena  taking  place  at 
the  Rayleigh  critical  angle,  Bertoni  and  Tamir  [2]  assumed  a^  ■ ag  be- 
cause the  pole-zero  pair  locations  did  not  vary  significantly  with  a^. 
Becker  and  Richardson  [4],  on  the  other  hand,  calculated  the  magnitude 
and  phase  of  R(kx)  for  a given  set  of  separate  loss  parameters,  assuming 
a constant  longitudinal -to- shear  attenuation  constant  ratio.  Finally, 


they  all  neglected  attenuation  in  the  liquid  medium. 

The  results  to  be  obtained  for  the  case  of  absorptive  media  in 
this  work  are  not  subject  to  the  above-mentioned  limitations.  Frequen- 
cy is  used  as  an  independent  parameter  and  therefore  leads  to  definite 
frequency-dependent  values  of  the  loss  parameters.  The  loss  parameters 
are  determined  from  attenuation  constants  experimentally  measured  for 
various  frequencies  as  described  in  Chapter  V.  This  procedure  was  se- 
lected in  order  to  allow  comparison  of  calculated  results  with  expert- 
mental  measurements  of  non-specular  reflectivity. 

C.  Absorptive  Reflection  Coefficient  near  Critical  Angies 

. 

After  incorporating  attenuation  into  the  reflection  coefficient, 

(one  can  now  determine  the  analytical  properties  of  the  absorptive  R (k  ). 

x 


In  the  following,  calculations  are  again  presented  for  water-stainless 
steel  and  water-Plexiglas  interfaces.  The  values  to  be  used  for  the 


r 

i 
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loss  parameters  are  tabulated  in  Table  2.  The  attenuation  constants 
were  measured  for  the  samples  which  were  used  in  experiments  described 
in  Chapter  V.  As  will  be  seen  below,  the  Introduction  of  attenuation 
leads  to  two  developments:  the  positions  of  the  singularities  are  dis- 
placed from  those  of  the  lossless  case;  and  the  magnitude  and  phase 
curves  are  different  from  the  lossless  ones,  especially  near  the  criti- 
cal angles.  At  these  critical  angles,  the  lossless  curves  (Figs.  3 and 
4)  are  not  smooth  but  the  absorptive  ones  become  smoothly-varying. 

The  first  development  is  Illustrated  by  tracing  the  loci  of  the 
singular  points  in  the  complex  k^-piane.  The  branch  points  ahd  their 
associated  branch  cuts  are  displaced  vertically  off  the  real  axis;  this 
displacement  increases  with  increasing  frequency-dependent  attenuation. 
Fig.  2.  Accordingly,  the  pole  and  zero  positions  also  change  with  the 
frequency.  Their  new  locations  are  determined  by  a computer  program 
developed  to  search  for  the  zeroes  of  the  numerator  and  denominator  of 
the  absorptive  Rfk^)  utilizing  the  Newton-Raphson  algorithm  (Appendix  B) . 
Figure  5 depicts  such  loci  of  the  absorptive  pole-zero  pair  as  the  fre- 
quency varies  up  to  kO  MHz  for  a water-stainless  steel  interface,  it  is 
noted  that  in  the  absorptive  case  the  angle  of  Incidence  at  which  the 
leaky  Rayleigh  mode  is  generated  depends  on  the  sound  frequency  since 
the  real  part  of  kp  now  assumes  a slightly  different  value  from  the  loss- 
less case. 

The  other  development  is  the  change  in  the  shape  of  the  magnitude 
and  phase  curves.  Again  the  magnitude  and  phase  of  R(kx),  which  now 
contains  the  loss  parameters,  are  calculated  for  the  water-Plexiglas 
and  water-stainless  steel  combinations.  In  the  water-Plexiglas  case. 


TABLE  2 


ATTENUATION  DATA 


f 

(KHz) 

®d 

(Neper/m) 

o 

s 

(Neper/m) 

*d 

as 

Stainless 

Steel 

1 

.027 

.025 

. OOOOA 

.000037 

2 

.229 

.330 

.00017 

. 0002A 

3 

.795 

1.A9 

.00039 

.00073 

6 

6.68 

19.6 

.00165 

.OOA83 

10 

32.1 

131 

.00A7A 

.019A 

15 

111 

590 

.0111 

.0583 

20 

269 

1720 

• 0199 

.1273 

25 

535 

39A3 

.0317 

.2333 

30 

936 

776A 

.0A62 

.3830 

35 

1503 

13770 

.0635 

.5823  ’ 

AO 

2265 

22620 

.0838 

.8369 

Plexiglas 

1 

7 

.010A 

7.2 

.0107 

2 

36.6 

.0271 

67.0 

.OA96 

3 

88. A 

• 0A36 

2A7 

.1219 

6 

AI6 

.1027 

2298 

.5670 

Water 

.00026* 

^The  attenuation  constant  for  water  ts  assumed  to  be  a linear 
function  of  frequency  and  therefore  Its  loss  parameter  is  considered 
to  be  a constant  for  the  frequency  range  used. 


the  difference  between  the  lossless  and  absorptive  Rfk^)  is  most  signi- 
ficant near  the  longitudinal  critical  angle.  Figure  6 shows  the  magni- 
tude and  phase  as  functions  of  the  incident  angle  in  the  vicinity  of 
the  longitudinal  critical  angle.  It  is  seen  that  the  magnitude  of  the 
reflection  maximum  drops  considerably  with  increasing  frequency  in  the 
magnitude  curve,  and  that  the  non-smoothly  varying  behavior  associated 
with  the  longitudinal  critical  angle  disappears  in  both  magnitude  and 
phase  curves.  It  is  also  observed  that  the  reflection  amplitude  is 
maximum  at  slightly  different  angles  of  incidence.  For  the  water-stain- 
less steel  interface,  the  change  in  these  curves  is  most  noticeable  in 
the  angular  range  containing  both  the  shear  and  Rayleigh  critical 
angles.  Figure  7 shows  the  familiar  "least  reflection  frequency"  ef- 
fect, which  has  been  reported  by  other  authors  [2,43,  arising  out  of  the 
inclusion  of  attenuation  into  the  lossless  reflection  coefficient.  This 
effect  is  characterized  by  the  appearance  of  a sharp  minimum  in  the  re- 
flection magnitude  curve  exactly  at  the  Rayleigh  critical  angle.  In 
Fig.  7.  this  sharp  dip  is  seen  to  vary  with  the  frequency  and  to  reach  the 
lowest  value  at  about  15  MHz.  It  can  be  shown  that  the  frequency  at 
which  the  reflection  amplitude  is  zero  at  the  Rayleigh  critical  angle 
corresponds  to  the  point  where  the  locus  of  the  zero  of  the  absorptive 
R(kx)  cuts  the  real  axis  in  the  complex  k^-plane.  It  should  be  pointed 
out  that  it  is  no  longer  tenable  to  maintain  the  long-standing  associa- 
tion of  the  longitudinal  critical  angle,  defined  by  Snell's  law,  with  a 
reflection  maximum  as  a result  of  the  inclusion  of  attenuation.  It  is 
also  true  that  the  shear  critical  angle  of  an  absorptive  liquid-solid 
interface  loses  Its  "critical"  character  primarily  because  it  is  too 


PHASE  CHANGE  (OEGREES)  REFLECTION  AMPLITUDE 


ANGLE  OF  INCIDENCE  (DEGREES) 

Figure  6.  Magnitude  and  phase  of  the  absorptive  R(k^)  near  the  longt 
tudinal  critical  angle  for  a water-Plexiglas  interface. 
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Figure  7.  Magnitude  of  the  absorptive  R(k^)  near  the  Rayleigh 
critical  angle  for  a water-stainless  steel  interface. 
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close  to  the  Rayleigh  critical  angle. 

In  this  chapter,  the  essential  features  of  the  plane  wave  nf lec- 
tion coefficient,  with  and  without  attenuation,  have  been  extensively 
examined  in  anticipation  of  its  determining  role  in  the  non-specular 
reflectivity  problem.  This  problem  itself  will  be  studied  In  the  re- 
maining chapters. 


CHAPTER  III 


THEORETICAL  FORMULATIONS  OF 
NON-SPECULAR  REFLECTIVITY 

This  chapter  presents  different  theoretical  formulations  designed 
to  explain  the  experimentally  observed  non-specular  reflection  phenomena 
which  arise  in  practice  when  a bounded  ultrasonic  beam  is  incident  onto 
an  interface.  All  formulations  to  be  discussed  here  are  constructed 
from  the  principle  of  spectral  representation  of  a finite  beam  but,  in 
general,  their  success  originates  largely  in  their  particular  method  of 
treating  the  complex  plane  wave  reflection  coefficient  which  plays  a 
central  part  in  these  formulations.  In  the  following,  the  common  me- 
thod of  representing  a bounded  beam  is  briefly  reviewed  before  the  var- 
ious formulations  are  discussed  together  with  methods  of  treating  the 
reflection  coefficient,  the  chapter  is  to  be  concluded  with  a critical 
evaluation  of  these  formulations. 

A.  Spectral  Representation  of  a Bounded  Ultrasonic  Beam 

Consider  an  Incident  beam  of  which  the  cross-section  In  the  (x,z) 
plane  Is  bounded  and  which  Is  uniform  along  the  y-direction,  that  is, 
there  is  no  spatial  y dependence.  With  respect  to  the  liquid-solid 
interface,  the  coordinate  system  Is  chosen  as  shown  in  Fig.  8,  where 
the  incident  beam  is  taken  to  be  Incident  at  an  angle  flj  and  Its  pro- 
file is  contained  by  an  effective  width  2w,  projected  onto  the  Inter- 
face as  2wq  given  by  wQ  • w sec  8|. 
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Figure  8.  A schematic  diagram  of  a bounded  ultrasonic 
beam  Incident  on  a liquid-solid  Interface. 
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It  Is  shown  [8]  by  using  Fourier  analysts  that  the  Incident  field  of  a 
beam  can  be  uniquely  determined  at  any  point  (x,z)  if  the  field  distri- 
bution is  known  in  any  plane.  Describing  the  acoustic  field  in  terms  of 

\ 

the  z-component  of  the  particle  displacement,  the  incident  field  can  be 
expressed  as  a sum  of  infinite  plane  waves  by  means  of  the  Fourier  in- 
tegral-transform pair 

-l4** 

UInc(x,z)  - (2u)  / V(kx)  expplk  x + k z)]  dkx  (5) 

+® 

V(kx)  - / U.nc(x,0)  exp(-lkxx)  dx  , (6) 

where  the  time  dependence  exp  (-ioit)  is  suppressed,  kx  is  the  x-compo- 

nent  of  the  wavenumber  vector,  V(kx)  is  interpreted  as  the  amplitude  of 

each  of  the  constituent  infinite  plane  waves,  and  the  z-ccmponent  of 

2 2 k 

the  wavenumber  vector  k2  Is  defined  by  kz  » (k  -kx)  . 

I.t  has  been  assumed  that  the  effective  beam  width  is  2w  and  the 

beam  strikes  the  interface  over  a width  2w  . For  a well-defined  beani, 

o 

neglecting  beam  spreading,  one  should  expect  contributions  to  the  inte- 
gral In  (5)  only  from  those  infinite  plane  waves  which  are  defined  by 

kj  - */wo  t kj  ♦ »/wo  (7) 

and  the  limits  of  integration  in  (6)  should  be 


-w  < x i ♦ w , (8a) 

o o 

where 
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kj  » | k | slnOj  ■ (w/v)  sin  0.  - (2»A)  sin  8j  . (3b) 

In  view  of  the  fact  that  the  integral  in  (5)  is  carried  out  along 
the  real  axis  and  that  is  varied  in  a narrow  range  around  k^ . it  can 

be  seen  from  (8b)  that  the  above  mathematical  representat ion  of  a bound- 
ed beam  entails  two  possible  physical  interpretations:  either  the 
bounded  beam  is  composed  of  an  infinite  number  of  plane  waves  having 
the  same  wavelength  but  incident  at  different  angles;  or  the  beam  con- 
sists of  parallel  plane  wave  rays  with  a small  perturbation  in  wave- 
length due  to  the  non-uniformity  of  the  sound  frequency.  In  practice, 
one  may  expect  both  physical  situations  when  the  experimental  procedure 
of  generating  the  sound  beam  is  brought  into  consideration.  In  this 
and  subsequent  chapters,  k^  is  treated  as  an  ultimate  independent  var- 
iable without  taking  into  account  the  physical  reasons  behind  its  var- 
iation. 

Extending  the  spectral  representation  of  a bounded  beam  to  the 
reflected  sound  beam,  one  can  also  represent  the  reflected  field  as  a 
superposition  of  infinite  plane  waves  [2].  The  reflected  beam  profile 
can  now  be  determined  by  integrating  over  individual  plane  waves  whose 
magnitude  is  multiplied  by  the  reflection  coefficient.  Thus  the  z-com- 
ponent  of  the  reflected  field  is  written  as 

U f (x,z)  - (2*)  ' / R(kx)V(kx)exp  [i  (kxx  + k^)]  dkR  , (9) 

where  R(kx)  Is  given  by  (1). 

At  the  interface,  z ■ 0,  (9)  reduces  to 
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uref (x»°)  - (2*f 1 / R(kx)V(k  )exp(!kxx)  dkx  . (10) 

According  to  (9),  the  reflected  beam  profile  Is  constructed  by 
the  interference  of  an  infinite  number  of  reflected  plane  waves,  which 
have  an  amplitude  V(kx)R(kx)  and  undergo  a phase  shift  introduced  by 
RCk^)  upon  reflection.  As  a result  of  this  very  method  of  constructing 
the  reflected  field  distribution,  deviation  from  the  incident  beam  pro- 
file is  expected  to  be  substantial  where  either  the  magnitude  or  the 
phase  of  R(kx)  changes  appreciably.  This  is  the  cause  for  the  earlier 
anticipation  that  non-specular  reflection  phenomena  would  be  most  visi- 
ble near  certain  critical  angles. 

Due  to  the  mathematically  complicated  form  of  R(kx)  in  (9),  the 
integral  In  most  of  the  cases  cannot  be  evaluated  in  closed  form.  In 
order  to  overcome  this  difficulty,  several  authors  have  developed  var-> 
ious  analytical  methods  of  simplifying  the  integrand  in  order  to  eval- 
uate this  integral.  The  remaining  sections  of  this  chapter  are  intended 
to  review  these  formulations  and  to  introduce  a numerical  integration 
method  together  with  a comparative  analysis  of  their  relative  usefulness. 
B.  Schoch  Treatment  of  Reflected  Beam  Profile 

On  the  basis  of  the  spectral  representation  principle,  Schoch  [l] 
sought  a theoretical  explanation  for  the  non-specular  reflection  phe- 
nomena. His  formulation  is  intended  to  investigate  the  reflected  beam 
profile  In  the  neighborhood  of  those  incident  angles  that  lead  to  total 
reflection.  Since  this  formulation  is  constructed  with  the  lossless 
reflection  coefficient  in  the  integrand  of  (9),  the  amplitude  of  R(kx) 
is  assumed  to  be  constant  for  the  considered  angles  of  incidence  and 
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IU  Phase  I,  assumed  to  vary  slowly  In  such  circumstances.  Under  these 
assumed  conditions,  the  plane  wave  reflection  coefficient  is  simplified 
by  a first-order  Taylor  series  expansion  to  become  [8] 

*(kxJ  » |R(k,)|  exp  [|(kx  - kt)*'(k,)]  , (n) 

»*here  #'  is  the  first  derivative  of  the  phase  *(k  ) of  R(k  ) and  both 

^ X 

the  magnitude  and  are  evaluated  at  k,.  Substituting  (1 l)  into  (9), 
one  obtains 

Uref(*.*)  « (2w)*,|R(k,)|«xp  C-Ik,4* (kf)] 

' Jm  V(kx}  »kx[*  ♦ ♦ ' (kj )J  + Ik 2z  } dkx  . (12) 


Comparing  (12)  and  (5),  It  is  found  that 

»l«(k,)|,xp  [-lk,*'(k,)]  UIne  [x  + ♦•(k,)j  . (13) 

Th*  expression  (13)  Indlcetes  that  the  beam  is  displaced  along 
the  interface  by  a distance 

AS  " “ • (Mi) 

for  the  reflected  beam  profile,  Schoch's  derivation  thus  predicts 
a lateral  shift  that  depends  solely  on  the  first  derivative  of  the  rela- 
tive phase  function  of  the  plane  wave  reflection  coefficient.  Schoch's 
result  has  been  shown  tn  many  cases  to  be  Incompatible  with  experimental 
observations  near  the  critical  angles.  Its  failure  can  be  traced  back 
to  the  facts  that  (l)  his  derivation  is  based  on  simplistic  assumptions 


Implied  in  the  behavior  of  the  lossless  reflection  coefficient,  and  (ii) 
it  does  not  consider  the  existence  of  singularities  of  R(kx),  that  is, 


the  singular  behavior  of  R(k^)  near  the  critical  angles.  In  retrospect, 
Schoch's  analysis,  in  its  simple  form,  has  shown  that  there  always 
exists,  although  insignificant,  a certain  amount  of  non-specular  reflec- 
tivity in  the  form  of  a beam  displacement  for  those  angles  of  incidence 
at  which  R(kx)  satisfies  the  assumed  conditions  on  magnitude  and  phase. 
C.  Bertoni  and  Tamir  Formulation  for  Incidence 
at  the  Rayleigh  Critical  Angie 

In  another  way  to  overcome  the  analytical  difficulty  caused  by  the 
complicated  form  of  R(kx),  Bertoni  and  Tamir  [2], taking  advantage  of  the 
existence  of  the  pole-zero  pair  associated  with  the  Rayleigh  critical 
angle  in  order  to  simplify  R(kx),  were  able  to  carry  out  the  Integration 
in  (9)  in  closed  form.  The  reflection  coefficient  is  expanded  in  terms 
of  a Laurent  series  with  only  the  most  significant  term  being  retained, 
giving  rise  to 

R<kx>  - F-nr  • (,5> 

x P 

A convenient  form  of  the  Incident  beam  is  a Gaussian  distribution 
of  half-width  w.  Hence  the  incident  sound  field  can  be  characterized  by 

Ujnc^x»°)  * MP  [“(x/w0)2  + lkjx]  » 06) 

which  Is  normalized  to  unity.  To  determine  the  reflected  beam  profile, 

one  needs  to  find  the  Fourier  transform  of  U.  (x,0)  in  the  k -plane. 

• nc  x 

Inserting  (16)  into  (6),  one  obtains 
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V(kx)  - A<o  exp  [-(kx  - kj)2(wo/2)2J  . 


07) 


Substituting  (15)  and  (17)  Into  (10),  the  reflected  field  dlstrt- 
button  at  z ■ 0 Is  found  to  be  [2] 


“reff*-0'  -CTT  Ulrc(x,0) 
• P 


(18) 


♦ kP-!-k°  [l  + wo(kp-kJ)exp(y2)erfc^)]  uinc(x»°). 

pi 


where  erfc(y)  is  the  complementary  error  function  of  argument 


w lm(k  ) x Os  * 8e (K  ) Jw 

t . °.  f.  - — ♦ i±J — -e-L a.  • 

O 


(19) 


For  incidence  at  the  Rayleigh  angle. 


kj  - Re(kp) 


and 


w lm(k  ) 
o P 


w_ 


the  incident  and  reflected  beam  intensity  can  be  described  by 


|U,nc^*°H2  " "P  [-2(x/wq)2] 


(20) 


and 


Iuref (x,0) |2  " IR1  ♦ R2R3|2  lujnc(x,°^  ^ » (2,) 

where 


w'V  lm(k  ) , 

Rj  ■ l E_  exp(yZ)  erfc(y)  . 


With  the  above  formulation,  BertonI  and  Tamir  have  provided  a 
satisfactory  model  in  describing  the  essential  features  of  the  non- 
specular  reflection  phenomena  for  incidence  at  the  Rayleigh  angle.  Ber- 
toni  and  Tamir's  model  demonstrates  that  the  reflected  beam  profile  ex- 
hibits a lateral  shift,  an  intensity  null  within  the  reflected  beam,  and 
a weakening  trailing  field  as  had  been  experimentally  observed [3].  A 
more  accurate  quantitative  description  can  be  obtained  by  this  mode)  if 
attenuation  in  the  media  is  included  in  the  formulation.  In  the  next 
chapter,  some  reflected  profiles  are  determined  numerically  from  Bertoni 
and  Tamir's  results,  including  attenuation,  for  comparison  with  others. 
The  computer  program  written  for  this  purpose  is  illustrated  in  Appen- 
dix C. 

By  considering  separately  the  two  terms  in  (id),  Bertoni  and  Tamir 


have  demonstrated  that  the  resulting  reflected  profile  was  constructed 
from  the  Interference  of  the  specularly  reflected  field  and  the  field  of 


leaky  Rayleigh  waves.  This  interpretation  is  quite  adequate  for  the  par- 


ticular case  of  non-specular  reflectivity  at  the  Rayleigh  angle,  but  is 
inappropriate  to  explain  the  same  phenomena  at  other  angles  of  incidence. 
This  also  suggests  some  limitations  of  the  Bertoni  and  Tamir  formula- 
tion: it  cannot  be  used  to  solve  the  non-specular  reflectivity  problem 
at  incident  angles  other  than  the  Rayleigh  angle,  and  it  is  not  applica- 
ble to  non-Gauss  I an  incident  profiles  in  genera). 

D.  Numerical  Integration  Method  of  Determining 
Reflected  Beam  Profile  for  All  Incident  Angles 

Oue  to  the  limitations  and  restrictions  of  the  Schoch  and  the  Ber- 
toni and  Tamir  models,  another  method  is  needed  to  investigate  the  non- 
specular  phenomena  at  all  angles  of  Incidence.  Using  the  spectral  re- 
presentation approach  for  a bounded  beam,  a numerical  integration  com- 
puter program  (Appendix  D)  is  developed  to  determine  the  reflected  beam 
profile  given  by  the  integral  expressed  in  (10).  The  program  is  built 
on  the  basis  of  the  Simpson's  integration  algorithm  applied  to  complex 
variables.  The  Integrand  is  rewritten  in  terms  of  dimenstonless  norma- 
lized variables  to  improve  accuracy  and  to  facilitate  the  input  of  physi- 
cal quantities. 

First, the  wavenumber  components  are  transformed  into  dimensionless 
quantities  through  the  following  new  variables 

\ * (v/'“,  k*  • 

A - 1 ♦ ia/2*  , 

Ad  - (v/vd)  (1  ♦ iad/2v)  , 

A#  - (v/vs)  (I  ■*.  las/2ir)  , 


Aj  - slnflj  . 
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The  square- root  terms  are  redefined  accordingly  as 
B - (A2  - A*)*  . 

Bd  - (Aj  - A2)*  , 

Bs  - • 

Then  the  plane  wave  reflection  coefficient  given  by  (1)  can  be 
rewritten  as 


A(AX) 


(2A«  ~ As>2  * 'Wl  - °AsVB 

(2A2  - A2)2  ♦ kWl  ♦ pA^/B  ' 


(22) 


In  order  to  allow  comparison  later  with  results  obtained  by  the 

Bertoni  and  Tamir  formulation,  the  incident  profile  used  is  chosen  to  be 

Gaussian  and  normalized  with  respect  to  unity  as  given  in  (16).  If  W 

denotes  the  projected  half  beam  width  measured  In  terms  of  the  sound 

wavelength,  i.e.  W » (w  /X),  and  X the  position  in  units  of  w , i.e. 

o o 

X «•  (x/w  ),  then  (10)  becomes 
o 


Uref(X) 


Al+Wk 

tt*)44*  / A(Ax)  V(Ax)  E(X,Ax)  dAx 


(23) 


where 

V(Ax)  - W exp  [-(itW)2(Ax  - A,)2]  , 


E(x,Ax)  - exp(i2irWAxX)  , 
and 


Wk  - 1/2W  . 
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The  program  Is  constructed  to  calculate  the  reflected  beam  pro- 
file In  the  form  of  (23)  derived  from  a normalized  Gaussian  incident 

profile.  A standard  accuracy  check  is  performed  on  the  iterative  inte- 

-8 

gration  process  and  it  is  found  that  an  optimized  accuracy  up  to  10 
can  be  attained  when  the  number  of  Intervals  reaches  60.  All  the  cal- 
culated results  shown  in  the  following  chapters  are  done  under  these 
conditions. 

E.  Comparative  Evaluation  of  Existing  Formulations 

In  this  section,  a brief  summary  is  given  of  the  procedures  used 
in  the  past  to  determine  the  reflected  beam  profile.  This  is  presented 
below  so  that  the  relative  advantages  and  disadvantages  of  the  existing 
methods  can  be  assessed. 

Under  fairly  strict  assumptions,  which  are  applicable  only  to  non- 
crltical  angular  ranges  that  insure  an  analytical  and  slowly-varying  be- 
havior for  the  lossless  R(k^),  Schoch's  simple  derivation  consists  of  an 
approximation  of  the  term  R(k^).  This  leads  to  a prediction  of  s beam 
displacement  without  the  need  to  carry  out  the  actual  integration. 

In  order  to  arrive  at  the  reflected  beam  profile  via  the  Bertoni 
and  Tamir  formulation,  the  following  steps  must  be  taken:  (i)  For  a 
liquid-solid  combination,  their  physical  properties  and  the  sound  fre- 
quency are  defined;  (ii)  The  Rayleigh  pole  and  zero  need  to  be  located; 

(ill)  Since  the  pole  and  zero  locations  change  with  frequency  when  the 
attenuation  in  materials  is  accounted  for,  their  loci  as  functions  of 
the  frequency-dependent  attenuation  In  the  complex  k^-plane  must  be 
found  if  a more  accurate  reflected  profile  is  desired;  (iv)  And  finally 
the  pole-zero  locations  are  inserted  into  the  closed  form  solution  given 

_ r 


- 37  - 


by  (21)  to  obtain  the  reflected  beam  profile  at  the  Rayleigh  critical 
angle. 

The  numerical  Integration  formulation  proceeds  in  a more  direct 
manner  without  going  through  the  mode-ident ifying  step  required  by  the 
Bertoni  and  Tamir  formulation,  ft  consists  of:  (i)  For  a give?!  liquid- 
solid  interface,  the  sound  velocities,  density  ratio  and  the  frequency 
together  with  its  associated  attenuation  parameters  are  determined; 

(ii)  These  values  are  then  considered  as  Inputs  to  a computer  program 
built  on  the  basis  of  the  Simpson's  integration  algorithm  providing  an 
accurate  profile  of  the  reflected  beam. 

It  is  obvious  that  the  contribution  of  the  Schoch  formulation  to 
the  theoretical  investigation  of  the  non-specular  reflection  phenomena 
is  the  most  limited.  It  fails  to  provide  a complete  physical  picture 
for  the  problem  near  the  critical  angles.  However,  it  should  be  realized 
that  Schoch 's  result  is  adequate  for  describing  non-specular  reflectivi- 
ty at  non-critical  angles  of  incidence. 

In  comparing  the  other  two  formulations,  one  sees  several  advan- 
tages of  the  numerical  integration  approach  over  the  Bertoni  and  Tamir 
analysis:  First,  the  numerical  integration  method  is  applicable  for 
all  angles  of  incidence  with  equal  ease  and  accuracy  and  thus  can  pre- 
sent a continuously  changing  profile  of  the  reflected  beam  as  a function 
of  the  Incident  angles;  Second,  one  does  not  need  to  go  through  the  time- 
consuming  process  of  identifying  the  Rayleigh  mode  in  the  case  of  Inves- 
tigating the  non-specular  reflectivity  effects  near  the  Rayleigh  angle; 
Third,  the  same  approach  can  be  readily  extended  to  study  more  subtle 
aspects  of  the  problem,  such  as  the  beam  spreading  effect,  non-Gaussian 
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incident  profiles,  and  the  internal  structure  of  the  bounded  beam;  And 
finally,  its  accuracy  is  sufficient  for  all  practical  purposes  and  can 
be  improved  with  better  possible  programming  techniques  unlike  the  other 
formulations  where  inaccuracy  is  Implicitly  built  in  by  various  approx- 
imations and  assumptions. 

In  conclusion,  it  should  be  pointed  out,  regarding  the  Bertoni 
and  Tamir  formulation,  that  thanks  to  Tts  ingenious  utilization  of  the 
Rayleigh  pole-zero  pair  It  is  able  to  present  the  non-specular  reflec- 
tion phenomena  at  theRayleigh  angle  from  another  useful  viewpoint:  an 
interaction  of  the  leaky  Rayleigh  field  and  the  specularly  reflected 
beam,  In  addition  to  the  process  of  interference  between  infinite  plane 
waves  as  implied  within  the  spectral  representation  framework. 


CHAPTER  tV 


REFLECTED  BEAM  PROFILE  CALCULATED 
BY  NUMERICAL  INTEGRATION  METHOD 

The  numerical  integration  method  and  other  theoretical  formula- 
tions dealing  with  the  non-specular  reflectivity  effects  were  presented 
and  evaluated  in  the  previous  chapter.  The  numerical  integration  method 
was  seen  to  possess  several  advantages  over  the  often  used  formulation 
developed  by  Bertoni  and  Tamir.  In  this  chapter,  these  advantages  are 
demonstrated  through  the  calculated  results  obtained  by  the  numerical 
integration  method  for  various  ranges  of  the  incident  angle,  some  of 
which  are  made  possible  only  by  this  numerical  approach.  The  principal 
parameters  affecting  the  reflected  beam  profile,  in  general,  are  the 
sound  frequency,  the  beam  width,  and  the  angle  of  incidence.  Since  the 
reflected  profile  is  found  to  have  distinct  patterns  for  different  ranges 
of  incident  angles,  the  non-specular  reflection  phenomena  are  to  be  in- 
vestigated separately  for  the  two  most  significant  ranges  near  the  Ray- 
leigh and  longitudinal  critical  angles.  The  shear  critical  angle  is  con- 
sidered as  part  of  the  Rayleigh  angular  range.  In  each  angular  range, 
the  reflected  beam  profile  is  determined  as  one  parameter  is  varied 
while  keeping  the  others  unchanged.  In  this  manner,  one  can  clearly 
discern  the  influence  exerted  on  the  profile  by  the  beam  width  and  the 
sound  frequency  through  attenuation  in  materials.  In  particular,  for 
the  case  of  incidence  at  the  Rayleigh  angle,  reflected  profiles 
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calculated from  the  Bertoni  and  Tamlr  formulation  are  compared  to  those 
obtained  by  the  numerical  integration  method. 

A.  Reflected  Profile  for  Incidence  near  Rayleigh  Critical  Angle 

tn  this  section,  the  reflected  profile  is  to  be  determined  for 
different  angles  of  incidence  near  the  Rayleigh  critical  angle  by  both 
methods,  numerical  integration  and  the  Bertoni  and  Tamir  method.  Since 
the  Bertoni  and  Tamir  method  is  an  approximate  formulation  providing  ac- 
curate results  only  for  the  Rayleigh  angle,  it  is  used  only  at  that 
angle  while  the  numerical  integration  approach  is  applied  for  all  angles 
of  interest.  The  profiles  produced  from  the  same  parametric  inputs, 
shown  in  Tables  1 and  2,  by  the  two  methods  are  compared  wherever  pos- 
sible. The  liquid-solid  combination  chosen  is  water-stainless  steel 
which  is  expected  to  exhibit  those  non-specular  reflectivity  effects 
near  the  Rayleigh  angle,  which  are  typical  for  many  liquid-metal  inter- 
faces. Similarly,  commonly  arising  values  of  the  sound  frequency  and 
the  beam  width  are  used  in  these  calculations. 

1 . Variation  of  Beam  Profile  with  Beam  Width 
A change  in  the  beam  width  is  expected  to  result  In  a change  in 
the  bounded  beam  internal  structure,  which  is  related  to  a change  in 
the  reflected  beam  profile  via  the  different  functional  behavior  of 
V(Ax)  and  E(A^)  in  the  integration  process  expressed  by  (23).  It  is 
realized  that  a good  understanding  of  these  interrelationships  would 
provide  a significant  insiqht  into  the  problem  of  constructing  a 
bounded  beam.  This  represents  another  major  task  beyond  the  scope  of 
the  work  undertaken  at  present.  Below  one  is  concerned  with  the  direct 
dependence  of  the  reflected  beam  profile  on  the  beam  width  itself. 
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The  reflected  beam  profile  ts  calculated  by  the  numerical  Integra- 
tion process  given  a fixed  frequency,  2 MHz,  exactly  at  the  Rayleigh 
angle  for  three  values  of  the  beam  width:  12  mm,  20  mm,  and  28  mm. 

The  results  are  plotted  In  Fig.  9 as  three  solid  curves  all  having  an 
Intensity  minimum  between  two  Intensity  peaks.  As  Indicated  by  these 
curves,  when  the  beam  width  becomes  larger,  the  complete  profile  grad- 
ually shifts  to  the  left  together  with  the  Intensity  null  while  the 
right  peak  gains  In  Intensity  and  the  left  one  loses  In  intensity.  The 
most  noticeable  change  comes  from  the  shift  of  the  intensity  null,  which 
Is  displaced  to  tne  left  by  17%  of  2wq  as  the  latter  Is  changed  from 
12  mm  to  20  mm.  For  Incidence  at  the  Rayleigh  angle,  the  same  reflected 
beam  profile  can  be  calculated  using  the  Bertoni  and  Tamir  formulation 
with  attenuation  taken  into  account.  These  results  are  presented  in 
the  same  figure  in  the  form  of  broken  curves.  Agreement  between  two  sets 
of  curves  Is  seen  to  be  reasonable  with  the  discrepancy  estimated  at  the 
intensity  peaks  to  be  about  2.5%  on  the  average.  The  outcome  indicates 
that  the  Bertoni  and  Tamir  formulation  Is  an  excellent  approximation  for 
the  description  of  the  non-specular  reflection  phenomena  at  the  Rayleigh 
angle  provided  that  the  physical  parameters  entered  are  sufficiently 
accurate.  The  numerical  integration  method  Is  more  dependable  in  view 
of  the  advantage  that  the  calculation  error  is  almost  non-existent. 

2.  Variation  of  Beam  Profile  with  Sound  Frequency 
The  sound  frequency  strongly  affects  the  reflected  beam  profile 
constructed  by  (23)  through  the  attenuation  parameters  in  R(A^)  and  the 
normalized  half  beam  width  W in  V(A^)  and  E(Ax).  The  relationship 
between  the  reflection  coefficient  and  attenuation  has  been  discussed  in 
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detail  in  Chapter  II.  The  numerical  integration  method  is  now  used  to 
calculate  the  reflected  profile  for  a fixed  beam  width  but  different 
sound  frequencies.  Figure  10  shows  these  profiles  for  a sound  frequen- 
cy of  1 MHz,  2 MHz,  and  6 MHz  and  a beam  width  of  20  mm  for  incidence 
at  the  Rayleigh  critical  angle. 

tt  is  observed  from  the  solid  curves  that  the  profiles  are  changed 
substantially  as  the  sound  frequency  changes.  Again  the  profiles  exhi- 
bit two  intensity  peaks  and  an  intensity  null.  As  the  frequency  is 
changed  from  1 MHz  to  2 MHz,  the  whole  profile  shifts  to  the  left,  with 
the  right  peak  becoming  higher  and  the  left  one  smaller,  while  the  null 
strip  moves  to  the  left  by  223;  of  2w^.  At  6 MHz,  one  sees  the  same 
changing  tendencies  except  that  the  left  peak  becomes  so  small  that  it 
almost  disappears.  In  the  same  figure,  the  reflected  beam  profiles  pro- 
duced from  the  Bertoni  and  Tamir  formulation  are  also  drawn.  Again  they 
are  seen  to  agree  with  the  results  of  the  numerical  integration  method 

and  the  approximation  error  involved  can  be  estimated  to  be  about  3«  at 

- 

the  intensity  peaks. 

3.  Variation  of  Beam  Profile  with  Angle  of  Incidence 

The  numerical  Integration  formulation  is  now  applied  to  invest!- 
gate  changes  in  the  reflected  beam  profile  as  the  angle  of  Incidence  is 
varied  about  the  Rayleigh  critical  angle.  This  is  the  case  where  the 
formulation  suggested  by  Bertoni  and  Tamir  is  not  fully  applicable.  It 
does  not  provide  accurate  results  for  angles  of  incidence  other  than  the 
Rayleigh  angle;  the  farther  from  it,  the  more  erroneous  the  Bertoni  and 


Tamir  results  become. 

Because  the  reflected  profile  exhibits  different  patterns  for 
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Figure  10.  Reflected  beam  profiles  for  a water-stainless  steel 
Interface  at  Rayleigh  angle  tnctdence  with  the  same 
beam  width  (20  mm)  but  different  frequency.  The 
solid  curves  are  calculated  by  the  numerical  Inte- 
gration method  while  the  broken  ones  are  produced 
by  the  Bertoni  and  Tamir  formulation. 
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different  frequency  ranges,  calculations  are  done  for  two  particular 
cases.  In  one  case,  the  frequency  and  beam  width  are  chosen  to  be  2 MHz 
and  20  mm;  in  the  other,  they  are  IS  MHz  and  20  mm.  In  both  situations, 
the  angle  of  incidence  is  varied  in  steps  of  0.25°. 

Figure  11  illustrates  the  2 MHz  case  with  five  reflected  profiles 
indicating  gradual  changes  as  the  angle  of  Incidence  is  slowly  increased 
above  the  Rayleigh  angle.  It  can  be  observed  from  the  figure  that  the 
changing  profile  is  characterized  by  the  following  developments:  (i) 

The  intensity  minimum  remains  essentially  at  the  same  position  but  its 
intensity  slowly  increases  until  it  is  completely  flattened  out  approx- 
imately at  1°  above  the  Rayleigh  angle;  (li)  The  right  peak  slowly 
loses  its  intensity  and  the  left  one  becomes  larger  until  they  gradually 
merge  into  one;  and  (iii)  Beyond  1°  above  the  Rayleigh  angle,  the  ori- 
ginal reflected  profile  turns  into  a Gaussian  one  whose  only  remaining 
non-specular  feature  is  a small  beam  displacement. 

The  results  for  15  MHz  are  presented  in  Figs.  12  and  13.  The  re- 
flected profile  now  is  found  to  have  a simpler  nature.  The  non-specular 
reflectivity  is  manifest  only  In  a lateral  shift  of  the  beam.  Figure  12 
shows  the  gradual  changes  in  the  reflected  beam  profile  as  the  angle  of 
incidence  becomes  less  than  the  Rayleigh  angle:  (I)  The  beam  displace- 
ment becomes  smaller;  and  (It)  The  maximum  peak  becomes  more  intense  as 
the  extent  of  non-specular  reflectivity  is  no  longer  noticeable  at  more 
than  2°  away  from  the  critical  angle.  The  beam  displacement  and  peak 
intensity  are  plotted  in  Fig.  13  as  functions  of  the  angle  of  incidence. 

In  both  situations,  additional  calculations  were  done  for  angles 
of  Incidence  smaller  than  the  Rayleigh  angle  in  the  2 MHz  case  and  larger 

I 
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Figure  H.  Reflected  beam  profiles  calculated  by  the  numerical 
Integration  method  for  a water-stainless  steel  Inter 
face  when  the  Incident  angle  Is  varied  slightly  from 
the  Rayleigh  angle  and  the  frequency  and  beam  width 
remain  the  same  (2  MHz  and  20  mm).  The  broken  curve 
Is  the  Incident  beam  profile. 
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Figure  12.  Reflected  beam  profiles  calculated  by  the  numerical 

Integration  method  for  a water-sta fnless  steel  Inter- 
face when  the  Incident  angle  Is  varied  slightly  from 
the  Rayleigh  angle  and  the  frequency  and  beam  width 
remain  the  same  (15  MHz  and  20  mm).  The  broken  curve 
is  the  incident  beam  profile. 
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ANGLE  OF  INCIDENCE  (DEGREES) 

Figure  13*  Beam  displacement  (A)  and  peak  intensity  (8)  of  the 
reflected  beam  profiles  near  the  Rayleigh  angle  for 
a water~sta inless  steel  tnterface  as  functions  of 
the  incident  angle.  The  sound  frequency  and  beam 
width  are  15  MHz  and  20  mm,  respectively. 
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tn  the  15  MHz  case.  The  results  are  found  to  be  almost  Identical  to  the 
profile  changes  described  above.  It  can  be  concluded  that  non-specular 
reflection  effects  near  the  Rayleigh  angle  are  most  pronounced  at  the 
Rayleigh  angle  and  become  less  so  as  the  angle  of  Incidence  departs  from 
the  critical  angle.  However,  the  effects  are  still  observable  as  far  as 
2°  away  from  the  Rayleigh  angle. 

B.  Reflected  Profile  for  Incidence  near  Longitudinal  Critical  Angle 

The  results  arising  out  of  the  investigation  of  the  absorptive 
plane  wave  reflection  coefficient  also  lead  to  the  belief  that  non- 
specular  reflectivity  may  exist  near  the  longitudinal  critical  angle. 

The  reflected  beam  profile  for  that  critical  angular  range  can  also  be 
determined  by  the  numerical  integration  method.  Again  the  reflected 
profile  is  calculated  for  different  parameters,  including  the  sound  fre- 
quency, the  beam  width,  and  the  angle  of  incidence. 

In  this  section,  all  calculations  are  done  for  a water-Plexiglas 
interface.  The  choice  of  this  liquid-solid  combination  is  made  in  ac- 
cordance with  earlier  observations,  discussed  in  Chapter  II,  about  the 
plane  wave  reflection  coefficient:  in  the  case  of  a water-Plexiglas 
interface,  R (k^)  exhibits  strong  variations  in  magnitude  and  phase  near 
the  longitudinal  critical  angle.  As  a result,  preliminary  calculations 
show  that  at  the  longitudinal  critical  angle,  the  reflected  beam  profile 
for  a water-stainless  steel  interface  is  only  slightly  displaced  by  about 

2%  of.  2w  , in  contrast  to  the  case  of  water-Plexiglas,  where  the  dis- 
o 

placement  becomes  much  greater  and  can  be  readily  observed  experimentally. 

1.  Variation  of  Beam  Profile  with  Beam  Width 


The  reflected  profile  is  calculated  by  the  numerical  integration 
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formulation  at  the  longitudinal  critical  angle  for  a fixed  frequency, 

2 MHz,  with  different  values  of  the  beam  width:  12,7  mm,  19,1  mm,  and 
25.4  mm.  The  non-specular  reflectivity  is  shown  by  a distinct  beam  dis- 
placement, which  is  calculated  to  be  14?,  13?,  and  11?  of  2wq  for 
12.7  mm,  19.1  mm,  and  25-**  mm,  respectively.  The  complete  profiles  for 
these  values  are  plotted  in  Fig.  14.  It  is  seen  that,  as  the  beam  width 
becomes  larger,  the  beam  displacement  decreases  and  the  peak  intensity 
inc reases. 

2.  Variation  of  Beam  Profile  with  Sound  Frequency 

The  reflected  beam  profile  is  calculated  for  different  values  of 
the  sound  frequency.  The  beam  width  is  constant  (19*1  mm)  and  the  angle 
of  incidence  is  the  longitudinal  critical  angle.  Figure  15  shows  the  re- 
flected beam  profiles  for  three  cases:  1 MHz,  2 MHz  and  6 MHz.  Again, 
the  only  non-specular  feature  seen  in  the  figure  is  a lateral  shift  of, 
the  reflected  beam:  11?,  13?,  arid  4!fc  of  2wq  for  1 MHz,  2 MHz,  and  6 MHz, 
respectively.  Because  of  lack  of  attenuation  data  for  other  frequencies, 
only  these  three  mentioned  profiles  are  calculated  and  thus  sufficient 
evidence  is  not  available  to  reach  a conclusion  regarding  variation  of 
the  beam  displacement  and  the  peak  intensity  as  functions  of  the  sound 
frequency.  However,  it  can  be  assumed  that  there  might  exist  a certain 
frequency  between  2 MHz  and  6 MHz  where  the  beam  displacement  and  inten- 
sity peak  would  reach  a maximum, 

3.  Variation  of  Beam  Profile  with  Angle  of  Incidence 


Figure  3 in  Chapter  U shows  that  near  the  longitudinal  critical 
angle  the  reflection  amplitude  is  strongly  dependent  on  the  angle  of 
Incidence,  increasing  from  40?  to  100?  reflection  within  less  than  2° 
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Figure  1 A . Reflected  beam  profiles  calculated  by  the  numerical 

integration  method  at  the  longitudinal  critical  angle 
for  a water-Plexiglas  interface  with  the  same  fre- 
quency (2  MHz)  but  different  beam  width.  The  broken 
curve  Is  the  Incident  beam  profile. 


Figure  15.  Reflected  beam  profiles  calculated  by  the  numerical  Inte- 
gration method  at  the  longitudinal  critical  angle  for  a 
water-PtexIglas  Interface  with  the  same  beam  width 
(19-1  mm)  but  different  frequency.  The  broken  curve  Is 


the  Incident  beam  profile. 
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change.  Expectedly,  the  reflected  profile  would  vary  appreciably  In 
this  angular  range.  The  numerical  Integration  method  Is  now  used  to  de- 
termine the  reflected  profile  of  a beam  of  2 MHz  frequency  and  19.  1 mm 
width  when  the  angle  of  Incidence  Is  varied  by  steps  of  0.5°  about  the 
longitudinal  critical  angle.  Some  of  these  profiles  are  shown  tn  Fig. 

16  for  incident  angles  between  0^  and  0^  + 3°.  it  can  be  seen  that  the 
variation  of  the  beam  displacement  and  peak  Intensity  are  substantial 
only  near  the  longitudinal  critical  angle. 

To  obtain  a more  detailed  knowledge  of  the  relationship  between 
the  beam  displacement  and  the  angle  of  incidence,  calculations  are  done 
for  smaller  increments  of  the  incident  angle.  The  calculated  results 
shown  by  the  solid  curves  in  Fig.  19  tn  Chapter  V indicate  that  there 
exists  a sharp  increase  in  both  the  beam  displacement  and  th*  peak  In- 
tensity approximately  at  0^  - 0.5°.  They  both  reach  a maximum  not  at 
0^  but  at  about  1°  above  it.  These  theoretical  predictions  are  consis- 
tent with  experimental  measurements  described  in  Chapter  V. 

It  Is  instructive  to  compare  the  above  profile  changes  to  those 
described  In  Fig.  12  for  the  water-stainless  steel  interface  at  15  MHz 
near  the  Rayleigh  critical  angle.  In  both  cases,  one  notes  that  the 
beam  displacement  increases  as  the  critical  angle  is  approached.  How- 
ever, the  peak  Intensity  rises  tn  the  water-? 1 ex !g las  case  and  decreases 
in  the  water-stainless  steel  case,  tn  view  of  these  observations  and  the 
known  behavior  of  R(k^),  one  may  arrive  at  the  following  association: 
the  peak  intensity  largely  expresses  the  variation  of  the  reflection 
amplitude  while  the  beam  displacement  changes  proportionally  with  the 
gradient  of  the  phase  shift  generated  upon  reflection. 


X(v*o/50) 

Figure  16.  Reflected  beam  profiles  calculated  by  the  numerical 
Integration  method  for  a water-Plexlglas  Interface 
v#hen  the  Incident  angle  Is  varied  near  the  longi- 
tudinal critical  angle.  The  sound  frequency  and 
beam  width  are  2 MHz  and  19. 1 mm,  respectively. 


CHAPTER  V 


EXPERIMENTAL  INVESTIGATION  OF 
NON-SPECULAR  REFLECTIVITY 

In  the  previous  chapters,  the  theoretical  aspects  of  the  problem 
of  non-specular  reflectivity  have  been  studied  in  detail,  first  by  ex- 
amining the  plane  wave  reflection  coefficient,  and  then  by  discussing 
various  formulations  for  determining  the  reflected  beam  profile.  Exper- 
imental evidence  of  the  non-specular  reflection  phenomena  has  been  re- 
ported by  several  authors  [3,7],  mostly  for  incidence  at  the  Rayleigh 
critical  angle.  With  the  numerical  integration  method  and  by  taking 
proper  account  of  attenuation  in  the  media,  the  non-specularly  reflect- 
ed profile  can  now  be  described  rather  accurately  for  ail  angles  of  in- 
cidence. 

A.  Scope  of  Investigation 

The  principal  objectives  of  the  present  experimental  investigation 
can  be  summarized  as  follows; 

(i)  Study  the  profile  changes  of  the  reflected  beam  when  the  angle 
of  incidence  is  slowly  varied  around  the  Rayleigh  critical  angle  for  a 
water-stainless  steel  interface. 

(il)  Monitor  the  changes  in  the  reflected  beam  profile  for  inci- 
dence near  the  longitudinal  critical  angle  onto  a water-Plexiglas  inter- 
face. A detailed  Investigation  of  the  beam  displacement  and  peak  Inten- 
sity is  undertaken  in  this  case  to  verify  earlier  theoretical  results. 
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B.  Experimental  Procedures 

The  development  of  the  experimental  procedures  used  here  is  cen- 
tered  around  two  standard  techniques:  the  Schlieren  imaging  technique 
to  visualize  the  path  of  a bounded  sound  beam,  and  the  pulse-echo  method 
to  measure  attenuation  in  materials.  The  experimental  work  proceeds  in 
the  following  steps: 

(I)  Determination  of  the  attenuation  parameters 

For  a specific  sample,  its  geometrical  dimensions  are  measured 
and,  together  with  its  mass,  are  used  to  determine  its  density  and  the 
length  of  the  travel  path  of  the  sound  pulse  through  the  sample.  Next, 
the  conventional  attenuation  constant,  a,  and  the  sound  velocity  are  de- 
termined by  the  pulse-echo  technique  for  a given  sound  frequency.  These 
results  are  used  to  calculate  attenuation  parameters  according  to  (1») . 
The  same  measurements  are  done  for  both  longitudinal  and  shear  waves. 

(II)  Study  of  the  profile  changes  for  a water-stainless  steel 
Interface  near  the  Rayleigh  critical  angle 

A Schlieren  imaging  experiment  is  set  up  to  observe  the  incident 
and  reflected  sound  beam  for  a stainless  steel  sample  immersed  in  water. 
The  steel  sample  is  25. A mm  thick  and  thus  can  be  considered  to  be  a 
half-space.  For  the  frequencies  used,  the  sample  surface  is  polished 
to  eliminate  effects  of  a rough  surface.  After  the  sound  frequency 
and  the  beam  width  are  chosen,  the  angle  of  incidence  is  slowly  varied 
about  the  Rayleigh  critical  angle.  The  continuously  changing  beam  pro- 
file is  first  recorded  by  a video  cassette  player  providing  a prelimi- 
nary survey  of  the  significant  features  of  non-specular  reflection  phe- 
nomena. The  angular  scanning  of  the  reflection  repeated  and  this  time 
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still  photographs  are  taken  of  the  complete  beam  profile  at  0.5°  incre- 
ments of  the  incident  angle,  including  the  critical  angle. 

(tit)  Study  of  the  profile  changes  for  a water-Plexig las 
interface  near  the  longitudinal  angle 
The  same  procedure  as  described  in  step  (ii)  is  employed.  The 
sample  thickness  again  is  25.4  mm  and  thus  would  qualify  as  a half-space. 
The  only  difference  is  that  the  angle  of  Incidence  now  is  varied  in  the 
neighborhood  of  the  longitudinal  critical  angle.  The  still  photographs 
of  the  beam  profile  are  also  taken  by  0.5°  steps.  In  addition,  a more 
detailed  verification  of  the  theoretical  results  predicted  by  the  numer- 
ical integration  method  is  attempted  through  a closer  measurement  of  the 
beam  displacement.  This  is  achieved  by  taking  more  photographs  of  the 
beam  profile  with  finer  increments  of  the  angle  of  incidence  and  measur- 
ing the  beam  displacement.  The  displacement  is  the  distance  between  the 
point  where  the  intensity  peak  of  the  incident  beam  strikes  the  inter- 
face and  the  point  of  reflection  of  the  intensity  peak  of  the  reflected 
beam. 

C.  Results  and  Analyses 
1 . Attenuation  Data 

Attenuation  data  for  Plexiglas  are  obtained  in  the  manner  describ- 
ed In  step  (i)  above.  The  sound  frequencies  used  for  both  longitudinal 
and  shear  measurements  are  1 MHz,  2 MHz,  and  5 MHz.  Interpolation  of 
the  measured  values  of  the  longitudinal  and  shear  attenuation  constants 
provides  the  complete  set  of  data  given  in  Table  2.  The  measured  sound 
velocities  in  Plexiglas  are  listed  in  Table  I.  The  interpolation  pro- 
cess is  achieved  by  fitting  the  experimentally  obtained  values  to  the 


curve  represented  by 


In  a ■ a Inf  - b , (24) 

where  f is  measured  in  MHz  and  a in  neper  per  meter.  For  the  shear  at- 
tenuation, the  coefficients  are  estimated  to  be  a * 3.218  and  b - 2.631, 
and  for  the  longitudinal  attenuation  a - 2.280  and  b ■ 2.659. 

Attenuation  data  for  stainless  steel  are  derived  from  measurements 
by  Fitch  [lo]  for  the  longitudinal  and  shear  attenuation  constants. 

These  values  are  again  fitted  by  the  same  method  and  are  extrapolated 
to  40  MHz.  The  coefficients  in  this  case  are  found  to  be  a * 1.614, 
b ■ 2.661  for  the  shear  attenuation  and  a * 1.334,  b - 2.626  for  the 
longitudinal  attenuation.  For  water,  the  attenuation  constant  is  taken 
to  be  a linear  function  of  the  sound  frequency  at  room  temperature  in 
the  frequency  range  of  interest.  As  a result,  the  attenuation  parameter 
of  water  is  frequency- independent,  as  shown  in  Table  2. 

2.  Profile  Changes  for  Water-Stainless  Steel  Interface 
The  reflected  beam  profile  is  observed  for  the  water-stainless 
steel  interface  near  the  Rayleigh  angle,  following  the  procedure  of 
step  (ii).  The  beam  width  and  the  sound  frequency  used  are  20  mm  and 
2 MHz,  respectively.  Figure  17  shows  five  Schlieren  photographs  of  the 
sound  beam  being  reflected  at  the  interface.  They  indicate  a reasonable 
agreement  with  the  theoretically  calculated  profiles  described  in 
Fig.  11.  As  the  angle  of  incidence  Is  Increased  from  1°  below  the  Ray- 
leigh angle  to  1°  above  it,  the  dark  strip  representing  the  intensity 
null  is  first  dimly  visible,  becomes  clearest  approximately  at  the 
Rayleigh  angle,  and  slowly  disappears  again  beyond  the  Rayleigh  angle. 


Figure  17.  Schlieren  photographs  showing  gradual  change  in  the  re- 
flected beam  profile  for  a water-stainless  steel  inter- 
face with  the  angle  of  incidence  being  (A)  - 1°, 

(B)  ep  - .5°,  (C)  (D)  er  ♦ .5°,  and  (E)  er  + i°. 


3.  Profile  Changes  for  Water-Plexig las  Interface 

For  the  water-Plexiglas  Interface,  the  reflected  profile  Is  ob- 
served as  the  incident  angle  is  varied  about  the  longitudinal  critical 
angle.  Non-specular  reflection  effects  are  evident  as  changes  in  the 
intensity  peak  and  the  beam  displacement.  The  sound  frequency  is  2 HHz 
and  the  beam  width  is  19.1  mm.  In  Fig.  )8  are  shown  four  photographs  of 
the  beam  profile  taken  from  Incidence  at  9^  ->  2°  to  incidence  at  8^  + 1°. 
A sharp  increase  in  the  beam  displacement  and  peak  intensity  can  be  de- 
tected for  an  incident  angle  slightly  less  than  6^.  Viewing  the  contin- 
uous picture  of  the  changing  profile  on  the  video  screen  indicates  good 
agreement  with  the  theoretical  predictions  given  in  section  IVB-3.  Mea- 
surements of  the  beam  displacement  on  photographs  taken  at  smaller  in- 
crements of  the  angle  of  incidence  are  plotted  together  with  calculated 
results  in  Fig.  19.  The  error  bars  on  curve  (C)  Indicate  the  estimated 
errors  in  measuring  the  beam  displacement  from  the  photographs.  The  dis- 
crepancy between  the  theoretical  beam  displacement  expressed  by  curve  (B) 
and  the  experimentally  measured  displacement  shown  by  curve  (C)  is  due 
to  the  fact  that  the  shear  wave  velocity  and  attenuation  are  difficult 
to  measure  to  a high  degree  of  accuracy.  However,  it  Is  significant  to 
note  that  both  calculated  and  experimental  curves  exhibit  the  same  var- 
iation pattern  for  the  beam  displacement  as  a function  of  the  angle  of 


incidence 
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beam  displacement  near  the  longitudinal  critical  angle.  The  angle  of  incidence 


The  theoretical  and  experimentally  measured  beam  displacements  are  denoted  by 


CHAPTER  VI 


CONCLUSIONS 

The  previous  chapters  have  been  devoted  to  a theoretical  analysis 
and  an  experimental  investigation  of  the  profile  of  a bounded  ultrasonic 
beam  reflected  from  a liquid-solid  Interface.  A numerical  integration 
method  has  been  employed  to  calculate  the  reflected  beam  profile  for  the 
complete  range  of  incident  angles.  Detailed  results  were  obtained  by 
this  method  with  different  parametric  conditions  for  two  particular 
cases:  incidence  near  the  Rayleigh  angle  onto  a water-stainless  steel 
Interface  and  incidence  near  the  longitudinal  critical  angle  onto  a 
water-Plexiglas  interface.  It  was  found  that  calculated  profiles  agreed 
reasonably  well  with  experimental  results  provided  by  the  Schlieren 
visualization  technique. 

The  investigation  has  led  to  the  following  major  results,  possible 
further  investigations  and  remaining  unsolved  questions: 

(i)  With  a rigorous  account  of  attenuation  in  the  media,  the  nu- 
merical integration  method  presents  an  accurate  theoretical  model  to  de- 
termine the  profile  of  a bounded  ultrasonic  beam  reflected  from  a liquid- 
solid  interface  as  a function  of  the  sound  frequency,  the  beam  width, 
and  the  angle  of  incidence.  The  most  significant  contribution  of  this 
method  is  that  it  allows  computation  of  the  profile  for  any  angle  of 
incidence  instead  of  being  restricted  to  non-specular  reflection  phe- 
nomena at  the  Rayleigh  angle.  This  aspect  is  demonstrated  by  the  result 
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obtained  in  the  present  work  for  the  case  of  a water-Plexiglas  interface 
in  which  non-specular  reflectivity  is  shown  to  exist  near  the  longitudi- 
nal critical  angle  and  can  be  described  in  quantitatively  satisfactory 
detai I . 

(ii)  Methodologically,  the  numerical  integration  method  affirms 
that  spectral  representation  Is  a sound  mechanism  to  express  the  bound- 
edness of  an  ultrasonic  beam.  It  entails  a physical  interpretation  of 
the  ben;'1  profile  as  the  interference  of  an  infinite  number  of  infinite 
plane  w.ii-es,  of  which  non-specuiar  reflectivity  at  the  Rayleigh  angle 
is  a particular  case  that  enables  one  to  view  the  reflected  profile  as 
a product  arising  out  of  a "resonance"  situation. 

(Ill)  In  view  of  the  practical  potential  of  the  numerical  inte- 
gration method,  the  present  research  can  be  extended  in  several  direc- 
tions. First,  the  numerical  approach  could  be  applied  to  study  the  non- 
specular  reflectivity  problem  for  other  structures  of  layered  media, 
such  as  a solid  plate  immersed  in  a liquid,  of  which  the  reflected  beam 
profile  has  been  established  for  plate-mode  incidence.  Second,  on  an 
equivalent  footing,  the  non-specular  transmittivi ty  phenomena  could  be 
investigated  in  a systematic  way  by  the  same  method  again  for  different 
layered  structures.  Third,  thanks  to  its  numerical  formulation,  the 
numerical  integration  method  could  be  modified  with  relative  ease  to  ac- 
commodate non-Gaussian  profiles,  which  may  be  encountered  more  often  in 
practice.  And  finally,  the  same  approach  could  be  extended  to  study  the 
so-called  beam  spreading  and  diffraction  effects:  "beam  spreading"  de- 
scribes a phenomenon  in  which  a bounded  beam  is  frequently  not  well- 
collimated  but  gradually  diverges  as  it  travels  in  a medium;  and 


•'diffraction"  indicates  a diffraction-like  beam  profile  when  observed 
away  from  the  active  surface. 

(Iv)  As  was  noted  earlier,  the  spectral  representation  principle 
Implies  that  a bounded  Incident  beam  is  constructed  from  a superposition 
of  Infinite  plane  waves  all  having  different  wave  vectors.  With  the  y 
component  being  eliminated  In  the  current  formulations,  the  difference 
between  the  wave  vectors  of  these  Individual  plane  waves  Is  manifest  In 
the  form  of  either  different  sound  frequencies  or  a perturbation  in  the 
angle  of  Incidence.  Reformulation  of  the  numerical  integration  method 
to  explicitly  consider  this  difference  may  help  understand  the  real  phy- 
sical structure  of  a bounded  ultrasonic  beam.  This  represents  a new 
problem  which  is  yet  to  be  solved  and  may  involve  some  aspects  of  non- 
1 inear  acoustics. 


APPENDICES 

FORTRAN  IV  COMPUTER  PROGRAMS  USED  IN  THEORETICAL  CALCULATIONS 
A.  Magnitude  and  Phase  of  Reflection  Coefficient 
jC  REFLECTION  COEFFICIENT  FOR  ABSORPTIVE  L/S  STRUCTURE 

1C  20  C WATER/STAINLPSS  STEEL 

C INPUT  PARAMEtERS 

P-3. 141592654 
V-1.48 
VO-5.84 
VS-3.13 

C CRITICAL  ANGLES  EVALUATED 

TANO-SQRT ( (V/VD)**2/( l-( V/VD)**2)) 

TANS -SORT I (V/V$)**2/(l™(V/VS)**2) I 

THETAD-ATANI TAND1 
THETAS-ATAN(TANS) 

. PS«-5*(VD**2-2*VS**2)/ I V0**2-V$**2) 

VR-VSM.  87+1. 12*PS ) / ll+PSI 
TANR-SQRTI ( V/VP ) **2/ ( l-( V/VR)**2) ) 

THETAR-ATANITANRI  _ „ _ 

TO-THETAD* 180/P 
TS-THETASM80/P 
TR-THETAR* 180/P 
WRITE  16,30) 

WRITE  16,35)  TD,TS,TR 
C INPUT  ATTENUATION  PARAMETERS 

05  READ  (5,40)  F,C2,C3 

IF  (C2.LT.0.0)  GO  TO  06 
' WRITE  (6,45) 

WRITE  (6,40)  F,C2,C3 

C SUBROUTINE  CALLED  TO  PRINT  OUT  ANPLITUOE  C PHASE 
I THETA-. 45  • _ 

STEP-I 

CALL  IMP( THETA, STEP,F,C2,C3I 
THET  A-THET  AD* 180/P 
' STEP-.Ol 

WRITE  (6,20) 

j CALL  IMP (THET A, STEP, F,C2«C3I 

THET  A-THET  AS*1 80/P 
WRITE  (6,25) 

CALL  I HP (THETA, STEP, F,C2,C3) 

THETA-THETAP*180/P 
i WRITE  (6,26) 

I..  _ CALL  I HP  (THET  A,  STEP,  F,C2,C3) _ 

! GO  TO  05 

| 06  STOP 

20  FORMAT!//  5K,*0ETAILE0  COMPUTATION  AT  LONG  CRIT  ANGLE* 

25  F0RMAT(//5X,*0ETAILE0  COMPUTATION  AT  SHEAR  CRIT  ANGLE* 

26  FORMAT! //5X, *OET AILED  COMPUTATION  AT  RAYLEIGH  CRIT  ANG 

30  F0RMAT(//5X,*THETA0*,9X,*THETAS*,SX,*THETAR*/) 

35  FORMAT (3E l 5. 4//) 


40  FORMAT  C10X,3FT0^6»  ' " 

45  FORMAT ( / 10X, • FREQUENCE*  *2X, *LONG  ATT* ,3X,* SHEAR  ATT*  I 


SUBROUTINE  IMP(THETA,$TEP,F,C2,C31 
COMPLEX  A(3),B(3),BP,BM,T1,T2«T3,R«CI 
1 _ OIMENSION  V(  3 1 1 C(  31  . . „ ..  _ . ... 

lC  INPUT  PARAMETERS 
< V<1)*1.48 

VC2)»5.84  _ ..... 

i VI3) *3. 13 

j P-3.141592654 

I . . DI-0. 12607 16  . . ......  .. 

CI-IO.O,l.O) 

Cl  I ) -0.00026/(2*  PI 

. CC  2I-C2/  ( 2*P  ) _ 

CI3I -C3/C2*PI 
DO  01  J«l,3 

AIJI-CVI I)/V(J))*CMPLX( l.O.CIJl)  

01  CONTINUE 

WRITE  16,151 

C EVALUATION  OF  REFLECTION  MAGN  .C  PHASE 

DO  02  K» 1,89 
21 -THETA-O.  45«-STEP*K 
AX-S IN( ZI*P/180I 
DO  03  L«lt3 
BP»AILI«-AX 

BM-AILI-AX  

PP-AT  A N2 (AIM AG IBP}, REAL (BP!)  . 

PM-ATAN2C AIMAGCBM) ,REAL I BN)  • 

IF  IPP.LT.O.O)  PP-PP«>2*P 
IF  ( PM.LT.O.O)  PN»PM*2*P 

Bl  L)-SQRTCCABSIBPI*CABS( BM||*CEXPfCI*f  PP*PMI/2 1 
03  CONTINUE  __ 

Tl-f AC3)**2-2*AX**21**2 
T2«4*BI2 l*B( 3)*AX**2 

! T3«0I*f A(3)**4l*B(2>/8(  1)  „ 

j R-ITUT2-T31/CTHT2+T3I 

RMAG-CABS(R) 

L . RPH-ATAN2(AtMAG<R!,REAL(R)l*l80/P ...  

IF  (RPH.LT.0.01  RPH»RPH«-360 
WRITE  16,10)  ZI,RMAG«RPH 

02  CONTINUE  _ 

RETURN 

10  F0RMATI3E15.4) 

15  FORMAT !6X, • INC  ANGLE* , 6X,* MAGNITUDE* ,10X, • PHASE*//) 
END  
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B.  Locations  of  Rayleigh  Pole-Zero  Pair 


DETERMINATION  OF  RAYLEIGH  POLE/ZERO  FOR  L/S  STRUCTURE 
20  C WATER/STAINLESS  STEEL 
*.  FREQ  RANGE  1-40  MHZ  - - — - 


COMPLEX  A(3),B(3I,CI *BPLS»BMNS 

COMPLEX  Dl, 02,03,001,002, 003, 0, DO, AP.APP  ‘ — * 

DIMENSION  VI  3) ,C 1 3) 

SOUND  VELS  AND  OTHER  CONSTANTS 

DATA  VI 1 1 , VI 2)  rv (3) /l. 48, 5 .84,3. 13/  : 

CI-10.0,  1.0) 

P-3.141592654 

DI-0. 12607 16"  ~ ‘ 

Cl  11*0.0 00 26 /1 2* P) 

WRITE  (6,40) 

'INPUT  ATTENUATION  PARAMETERS  "1 

READ  (5,35)  F,C2,C3 
C(2)-C2/(2*PI 

CI3J-C3/ I 2*P  I 

IF  (CI2I.LT.0.0)  GO  TO  31 

00  LOOPS  SCANNING  INITIAL  POINTS  START  HERE 

DO  20  1-1,6  ‘ 

00  15  J»l,6 

UP-0. 5095*1*0. 0005  • 

WP«- 0.  01  *J*  0.005  ' r 

AP*CMPLX(UP»WP) 

IC-l 

NORMALIZED  K • S~AND  KAPPAS  ' 

CONTINUE  • 

DO  10  K«l,3  I 

A(K)«(V( l)/V(KI l*CMPLX( 1.0«C(K)|  j 

BPLS-A (K I ♦ AP 
BMNS-A  IK  l-AP 

PPLS-ATAN2IAIMAGIBPLS) ,REAL(BPLSI) ' T 

PMNS-AT  AN2 I A I MAGI BMNS ) ,RE AL I BHNS ) ) 

IF  I PPLS.LT. 0.01  PPLS-PPLS*2*P 

IF  IPMNS.LT.0.01"  PMNS-PMNS*2*P 

BIK)«SQRT(CABS(6PLSI*CABStBMNS))*CEXP(C1*(PPLS+PMNSI/2] 

CONTINUE 

NUMERATOR  AND“ITS  DERIVATIVE * 

POLE  SEARCH  • 

Dl-IA(3I**2-2*AP**2)**2 

02- 4*B(2)*B(3)*AP**2 

03- 01*  81  21*1  A(  31**4)  /BID 
D-01*02*  D3 


1 
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. DABS«CAB$(D»  

DOl-8*APM2*AP**2-A(3l**2l 

002-8*AP*B(2l*B(3l-4*(AP**3)*(B(3l/B(2!+B(2!/B(3n 
DD3-(D1*AP*( A( 31**4! /(A(l)**2-AP**2) )*(B(2!/B( l I— BC 1178(21 1 
OD-OOl  ♦002003 
C NEWTON-PAPHSON  STEPS 

APP-AP— 0/00  ' 

OIFFR»REAL(APP-API 
_ 01FFI-AIHAGIAPP-AP1 

IF  (UBS(0IFFR)»ABS(0IFFI)!.I:e;1.E-07I  GO  TO  25 

AP-APP 

IC-!C*l 

IF  ( IC.GTiiSI"  GO*  TO  15  “ 

GO  TO  05 
15  CONTINUE 

20  CONTINUE  7 " 

GO  TO  30  * 

25 WRITE  (6,501  F,C2,C3,APP  < 

GO  TO  01 

30  WRITE  (6,551 

GO  TO  Oi  * 

31  STOP 

35  FORMAT  (10X,3FI0.6I 

60  FORMAT  (/4X,  »FREQ» *6X,  *tONG  ATT* ,3X, • SHEAR  ATT* ,4X, *REALPAI 
17X, * IMAGPART*/! 

50  F0RMAT(3E12«4,2E 15*71 

55  FORMAT  (/5X, *CANT  FIND  POLE*/!  * *1 

ENO 
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C.  Reflected  Profile  Calculated  from  Bertonl  and  Tamir  Results 


C 

C 

z~ 


REFLECTED  BEAR  FRCFIIE  FCR  ABSGRPT IVE  L/S  SYSTEM 
BEAR  WICTR  RANGES  FRCP  6 RR  TC  14  MR 

“AT  ~2~  ARC-*  2 

CCRPLEX  AP,A2,A1,R2 

RATER  VEL*KR/SEC,8EAM.  WICTP-MM.FREC-MhZ 

V-1746 — 

P-3«1415S2£54 

READ  IN  FPECePCLE  AND  ZERC 

“01  REACT  5720  J~ F • AFRYAP 1 • A ZR  VA2 1 

IF  (F.LT.C.C)  GC  IC  C6 
WRITE  (6,251 

WRITE  T672CI  F7APR,AP1,A2R7A21 

DC  06  J-1,5 
W-6.0«2.G«(  j-1* 

WRITE  (6*30 1 "W : 

WRITE  (6.10) 

AP-CNPLXI APR »AP  1 1 

A7-CRPCX(A7R7X7I J : 

Wl-W/SCRK  l-C-AFP**il 
ALFAW»P*F*API<k 1/V 

K0PMAL*17EC~V1NX~ANX~VREFL 

CO  05  1-1, 3CC 
X— 2*Wl*I*Wl/50 

G AM  H A- AL  F A W-X  / W 1 

Rl« (APR— A 2 )/(APR— API 
R2- ( AP— A2  J/IAP— APR1 

Rf-1— SCRTIP  }<ALFAW*EXP(GARRA**2)4ERFC(GARHA) 

VABS-CA8S(RWR2*R  > 

V1-EXP1-2MA7W1  J«»2I 

VR=YI*VA8S»*2 T 

WRITE  (6,15)  I, VI, VP 
CONTI RLE 

CCKTTNUF 

GC  TO  01 
STOP 

FORMAT— 172X7*  3riW075Cn  75X7’ V INC  •711X.*\REFL»/I 

FCR RAT  (110,2E15.6) 

FORMAT  ( 1CX ,5F 1C«7 I 

FORMAT-I /r3X7*FR  EG^76X7*  APR«77X7*AP  PTIXT*  AIR  • 77*7 
FORMAT  </10X,«W«*,F5 .1/1 
END 


H 


AZI 


i 
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D.  Reflected  Profile  Determined  by  Numerical  Integration  Method 


REFLECTED  BEAM  PROFILE  FOR  ABSORPTIVE  L/S  SYSTEM 
BY  OIRECT  INTEGRATION  USING  SIMPSON'S  RULE 
NEAR  THE  LONG  CRITICAL  ANGLE 

FOR  HATER/PLEXIGLASS  INTERFACE^  _ 

COMPLEX  Cl »FNE « FNR*  SUM2*  SUM4* VR  *VR1 • VR2  » VR3 
OEFINE  TWO  FUNCTIONS  IN  THE  INTEGRAND 

" -FNV( AX 1-EXPt— (CF*( AX— AI J/2 )**2I~~  

FNE I AX )=CEXP(CI*CF*AX*XNI 
INPUT  PARAMETERS 
' * N-60 
F-2.0E06 
V«1.‘8E03 
P-3.141592654 

CI-(0.0*1«0)  , 

THETA*32.69  . 

AI»SIN I THETA* P/180.0 1 
00  18  M-i,3 

W-2. 54*0.25*18*11/2.0  

W-W*0.01  

COMMON  FACTORS 
WX-W/SQRTI  l.C— A 1**2  ) 

CF-2*P*F*WX/V  ~ ~“ 

WK-P/CF 

WRITE  (6*20)  THETA, F,W 
"WRITE" 16*25)  N 
WRITE  (6,30) 

00  LOOP  TO  DETERMINE  UINC  ANO  UREF  AT  EACH  XN 

DO  "15  1-40*200 

X— 2*WX*I*WX/50 
XN-X/WX 

INTEGRATION  PARAMETERS  ANO  SUMMING  VARIABLES 
SUM4- 10*0*0.0)  , 

SUM 2- (0 .0*0.0 ) 

‘ FN-N : — 

H«2.0*WK/FN 

AX1-AI-WK 

AX2*  A I ♦ WK 

AX-AX1+H 

J«1 

SUMMIN6- ITERATIONS T 

SUM4-SUM4+FNR ( AX l*FNV( AX) *FNEI AX  I 
SUM2— SUM2*FNR( AX*H) *FNV( AX*H)*FNEf AX*H) 

IF  ( J.GE^N-3)  ~GG"T0"~1G : 

J-J42 
AX— AX*2*H 

GO  TO  05 


- 72  - 


^HfO  VR 1-FNR  C AX  1 » *FNV  I AX  1 ) *FNE  « AXl  I 
VR2*FNR(AX2l*FNVf AX2t*FNE(  AX2I 

VR3*FNR(AX2-H|*FNVCAX2-H|#FNE(AX2-H) 

VR-IH/3»M4*SU>K*2*SUM24VR1*4*VR3*VR2I 
UREF* ( CF**2/ ( 4*P I >*CABS ( VR I **2 

UINC-EXPI— 2*XN**2I 

WRITE  ( 6*35?  I.UINC.UREF 
15  CONTINUE 

~'18  CONTINUE  - 

STOP 

20  FORMAT  t/2X»  *THETA*« «E10.4» 2XV  «FREQ«*  *E8.2t2Xv  'WIDTH** 

25  ' FORMAT  ( /5X« * INTERVAL  NUMBER* *t 18)  ' ' 

. 3C  FORMAT  I/8X, • l • ,8X* • UINC* , 12X, « UREF « / ) 

35  FORMAT  (2X, I 10, 2E15.6 l 

END 


COMPLEX  FUNCTION  FNRIAX) 

COMPLEX  Cl«A(3)t8(3)*8P»BM*Tl«T2fT3 
DIMENSION  VI 3 1 * ATT (3) 

P—3.141592654 

01*0.8399832 

CI*IC.0,1.0) 

C SOUNO  VELS  ANO  ATTENUATION  CONSTANTS 

1K1M.48  

VI2I-2.74 

VC3I-1.35 

' ATTtU»0.00026/T2*pr 7 

AT TI2) *0.0271/1 2* PI 
ATT I 3 1*0.0496/1 2*P I 
-c  00  L0Qps  F(JR  Kt$  AND  KAPPAS 
00  50  K»i,3 

A(K)*IV(1)/VIKI l*CMPLXI l.Ct ATT (Kl I 
50  CONTINUE  * 

00  55  4.-1  *3 
8P»ACLI*AX 

8M-AILI-AX  * * ' ~ 

PP—ATAN2I AIMAG (8P) « REAL (BPI ) 

PM-ATAN2C AIMAG I BM I * REAL  IBM ) I 

“ IF  (PP.LT.O.OI  PP»PP+2*P  ' 

IF  (PM.LT.O.O)  PM*PP424P 

81 L l*SQRT(CABS( BPI*CABS I BMl l*CEXP ICl*( PP4PMI/2 I 

[ 55  CONTINUE 

C FUNCTION  EVALUATION  STEPS 

Tl-fAC3»**2-2*AX**2)**2 

“ T2*4*8(2)*8(3I*AX**2 

T3»0I*fA(3 14*41*8121/811) 
FNR-IT14T2-T3I/IT14T24T31 

RETURN 

ENO 
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